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NOTICE TO AUTHORS 


1. Communications.—Papers must be communicated to the Society by a Fellow. They 
should be accompanied by a summary at the beginning of the paper conveying briefly the 
content of the paper, and drawing attention to important new information and to the main 
conclusions. ‘The summary should be intelligible in itself, without reference to the paper, 
to a reader with some knowledge of the subject; it should not normally exceed 200 words 
in . Authors are requested to submit MSS. in duplicate. These should be 
eyped ws double spacing and leaving a of not less than one inch on the 
left-hand side. Corrections to the MSS. be made in the text and not in the 

Unless a paper reaches the Secretaries more than seven before a Council 
meeting it will not normally be considered at that meeting. By Council decision, MSS. of 
accepted papers are retained by the Society for one year after publication; unless their 
return is then requested by the author, they are destroyed. 


2. Presentation.—Authors are allowed considerable latitude, but they are requested to 
follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or in 
a numbered list at the end of the paper. Each reference should give the name and 
initials of the author cited, irrespective of the occurrence of the name in the text (some 
latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively :— 

A. Corlin, Zs. f. Astrophys., 15, 239, 1938. 
H. Jeffreys, Theory of Probability, 2nd edn., section 5.45, p. 258, Oxford, 1948. 


3. Notation.—For technical astronomical terms, authors should conform closely to the 
recommendations of Commission 3 of the International Astronomical Union (Trans. 
1.A.U., Vol. V1, p. 345, 1938). Council has decided to adopt the 1.A.U. 3-letter abbrevi- 
ations for constellations where contraction is desirable (Vol. IV, Ps 221, 1932). In general 
matters, authors should follow the recommendations in Symbols, Signs and Abbreviations 
(London: Royal Society, 1951) except where these conflict with I.A.U. practice. 


4. Diagrams.—These should be designed to appear upright on the page, drawn 
dea twice the size required in pa and prepared for direct photographic 


reproduction except for the lettering, ns should be inserted in pencil. 
Legends should be given in the pr awe tage cating where in the text the 
figure should appear. Blocks are retained by Society for for 10 years; unless the author 
requires them before the end of this period they are then destroyed. 


5. Tables —These should be arranged so that they can be printed upright on 
the page. 


6. Proofs.—Costs of alterations exceeding 5 per cent of composition must be borne by 
the author. Fellows are warned that such costs have risen sharply in recent years, and it 
is in their own and the Society’s interests to seek the maximum conciseness and simplifi- 
cation of symbols and equations consistent with clarity. 


7. Revised Manuscripts.—When papers are submitted in revised form it is especially 
requested that they be accompanied by the original MSS. 


Reading of Papers at Meetings 


8. When submitting papers authors are requested to indicate whether they will be 
willing and able to read the paper at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 


g. Postcards giving the programme of each meeting are issued some me before the 
meeting concerned. Fellows wishing to receive such cards whether for Ordinary 
Meetings or for the Geophysical Discussions or both should notify the Assistant Secretary. 
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MEETING OF 1954 NOVEMBER 12 
Dr J. Jackson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 

Frederick Benario, 151-10 State Street, Flushing 54, New York, U.S.A. 
(proposed by R. M. Baum); 

Mary Erwin Martin, 16 Rutland Gate, London, S.W.7 (proposed by 
F. J. M. Stratton); and 

Leslie Douglas Scott, County Laboratory, Dorchester, Dorset (proposed by 
A. F. O’D. Alexander). 


The election by the Council of the following Junior Members was duly 
confirmed :— 

Derek John Blundell, Imperial College of Science and Technology, London, 
S.W.7 (proposed by J. M. Bruckshaw); 

Eric Gray Forbes, B.Sc., Imperial College of Science and . Technology, 
London, S.W.7 (proposed by E. Finlay-Freundlich); 

Philip James Message, B.A., Gonville and Caius College, Cambridge 

_ (proposed by B. G. Tunmore); 

Richard Edward Small, University College, London, W.C.1 (proposed by 
C. W. Allen); 

David Edmund Smith, University College, Durham (proposed by F. 
Sheldrake); and 

Malcolm Walsh, University College, London, W.C.1 (proposed by C. W. 
Allen). 


Ninety-four presents were announced as having been received since the 
last meeting. 


MEETING OF 1954 DECEMBER 10 
Dr J. Jackson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 

Michael Abdulahad, Mesopotamia Observatory, Basrah, Iraq (proposed by 
A. K. Das); 

Christopher Rodney Armstrong, B.Sc., Guy’s Hospital, London, S.E.1 
(proposed by R. L. Waterfield); 
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Attia Abd-el-Salam Ashour, Faculty of Science, the University of Cairo, 
Egypt (proposed by V. C. A. Ferraro); 

Ivan Atanasijevi¢, The Physical Laboratory, Faculty of Sciences, Belgrade, 
Yugoslavia (proposed by M. Ryle); 

Sudhir Kumar Dutt, M.Sc., D.I.C., Imperial College of Science and 
Technology, London, S.W.7 (proposed by J. M. Bruckshaw); 

James Fogarty, 20 North Queen Street, Belfast, Northern Ireland (proposed 
by D. R. Bates); 

Kenneth W. Gatland, 431a Chertsey Road, Whitton, Middlesex (proposed 
by P. A. Moore); 

Morley Edwin Good, 18 Wardley Avenue, Little Hulton, Walkden, 
Manchester (proposed by E. Burgess); 

Charles West Graham, M.Inst.B.E., M.I.N., 3 Kathleen Avenue, Cleethorpes, 
Lincolnshire (proposed by F. E. Townend); 

Robert Hardie, Lowell Observatory, Flagstaff, Arizona, U.S.A. (proposed 
by S. Chandrasekhar) ; 

Leroy Freame Henderson, 70 Park Road, Middle Park, Melbourne, Australia 
(proposed by C. W. Allen); 

Theodore Edwin Houck, Washburn Observatory, University of Wisconsin, 
Madison, Wisconsin, U.S.A. (proposed by W. H. van den Bos); 

Henry George Hughes, M.Sc., F.Inst.P., M.I.E.E., 6 The Dale, Widley, 
Portsmouth, Hampshire (proposed by C. R. Burch); 

George W. Hutchinson, M.A., Ph.D., 12 Morven Road, Bearsden, Glasgow 
(proposed by T. R. Tannahill); 

Alan Maxwell, M.Sc., Ph.D., Jodrell Bank Experimental Station, Lower 
Withington, Macclesfield, Cheshire (proposed by A. C. B. Lovell); 

Sheila Constance Minion, B.Sc., 30 Fairfax Road, Derby (proposed by 
W. H. McCrea); 

Donald Mugglestone, University College of the West Indies, Jamaica, 
British West Indies (proposed by S. Chandrasekhar) ; 

William Nicholson, B.Sc., Royal Greenwich Observatory, Herstmonceux 
Castle, Sussex (proposed by D. H. Sadler); 

Jagdeo Singh, M.Sc., D.I.C., Imperial College of Science and Technology, 
London, S.W.7 (proposed by J. M. Bruckshaw); 

Henry Joseph Smith, Harvard College Observatory, Cambridge, Massa- 
chusetts, U.S.A. (proposed by R. H. Stoy); 

Peter Smits, B.Sc., 13 Cluny Road, Forest Town, Johannesburg, South 
Africa (proposed by A. D. Thackeray); 

Gladys Ethel Stone, 5 Belle Vue Road, Rodwell, Weymouth, Dorset 
(proposed by P. A. Moore); 

Antony Geoffrey Strutt, O.B.E., B.Sc., 49 Maintenance Unit, R.A.F., 
Colerne, Wiltshire (proposed by A. F. Collins); 

Jaakko Tuominen, Ph.D., Radio-astronomy Station, Siltavuorenpenger 20, 
Helsinki, Finland (proposed by F. J. M. Stratton); 

Chester B. Watts, U.S. Naval Observatory, Washington 25, D.C., U.S.A. 
(proposed by R. d’E. Atkinson); and 

Roy Bernard Worvill, M.Sc., Innisfree, The Leys, Chipping Norton, 
Oxfordshire (proposed by H. Wildey). 
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The election by the Council of the following Junior Members was duly 
confirmed :-— 


Charles David Austin, University College, London, W.C.1 (proposed by 
C. W. Allen); 

Michael Gadsden, Imperial College of Science and Technology, London, 
S.W.7 (proposed by F. J. Hargreaves); and 

Michael Guest, 122 Perrywood Road, Great Barr, Birmingham 22 (proposed 
by P. A. Moore). 


Ninety-seven presents were announced as having been received since the 
last meeting, including :— 


B. E. Markaryan, Atlas of known stellar clusters of varying types (presented by 
V. A. Ambartsumyan); 

A. W. J. Cousins, Standard magnitude sequences in the E regions (presented 
by the author); 

E. G. R. Taylor, The mathematical practitioners of Tudor and Stuart England 
(presented by the Cambridge University Press); and 

The Hon. Sir Charles A. Parsons—born 100 years ago on 13th Fune, 1954 
(presented by Messrs C. A. Parsons and Co., Ltd). 





SOLUTIONS OF THE NEGATIVE EMDEN POLYTROPES 


A. W. Rodgers and D. M. Myers 


(Communicated by the Commonwealth Astronomer) 
(Received 1954 September 7) 


Summary 


‘Tabular and graphical solutions of the negative Emden polytropes arising 
in the theory of spherical nebulae and clusters are found using a differential 
analyser. 

Models giving the light and density as a function of radius projected 
against the celestial sphere are constructed for five cases of the assumed 
mass-function. ‘The results indicate slight dependence of central intensity on 
this mass-function. 





1. Derivation of the equations.—A general theory of globular clusters has been 
presented in a recent paper by Woolley*, taking both dispersion in mass and trun- 
cation of velocity into account. The negative Emden polytropes arise as the 
equations governing the distribution of light and density with radius, and it is the 
purpose of this present paper to consider the effect on these distributions of 
arbitrary variations in the assumed mass-function. 


Woolley assumed a mass-function of the form C,, =const. x m*e~™ and then 
showed that the central mass-function &(m) is related to C,, by the equation 


Eq(m) a Cu Lm = erf, (% Vin)/ | Cr, Sm = erf, (% V Gm) dm, 


“z 
where erf,(x) is written for A | x2e-* dx. 
Valo 


Taking the velocity distribution to be fully Maxwellian, i.e. vg = 00, € (m) then 
becomes 


- oO 
—3/2 —3/2 
Contin 9] | Cyt 98 dm 
0 


Substituting for the assumed C,,, we have &)(m)am*%?e-°", The mean 


mass m is defined by m= | m2C_, dm | m 12 C, dm, =(a+ 4)/b and letting 
0 


0 


p=m m, 
C,,=const. x pt e-(@+H2) | 
Assuming, in the general case, equipartition of energy among the low-velocity 
stars, the density p reduces to the form 


const. x i; p2C,e“"’ du, where ¢ is the potential. 





*R. v. d. R. Woolley, M.N., 114, 191, 1954. 
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i) 
Then »=p/p)=const. x (a+$+y)-@t)® J x*-l2 e-*dx and since »=1 when 
0 


ys =o at the centre of the cluster, 


a+ a+1/2 
= (Sy) 
op 


If y=1+ —.,, then »=y “+, which, together with Poisson’s equation 
¥y ree yy" g q 


ads 
S ae =—4nrlI'p 
+) ’ 


gives " - (= =) =n, 
where z=r x (87Ip)8m)'?. 
Thus, if x=(a+ })"?z, we find that 


d dy 
eae 2s) me yay (2 +12 
z( 2) dase (1) 


with the boundary conditions y= 1, 2 =0 at x=0. 
This equation is called the negative Emden polytrope of index n =a+ 4 and has 
been solved numerically for the case a= 2-5 in Woolley’s paper. 

In this paper the index m was taken to have the five values, 2, 2-5, 3, 3-5 and 4, 
and the solutions of equation (1) were found using a differential analyser* con- 
structed by the Mathematical Instruments Section of the (Australian) Council for 
Scientific and Industrial Research Organization and housed in the Department 
of Electrical Engineering, University of Sydney. 

Letting y=rx', the asymptotic solutions of these equations are obtained by 
2(" + 3) 


substitutingin(1), givingr =p-/0+™, p = (n+1) 


andt=2/(1+m). Theasymptotic 


form of the density is then 
p=V M == pMl(L+n) y-Bn/(1+n)_ (2) 


In some casesit may be more convenient to compare the light distribution rather 
than the density with observation, and in order to obtain the light distribution we 
need to know the mass-luminosity law for stars in, say, globular clusters. 

We will assume that the law is L « M*°, though it is open to question whether 
this is applicable to Baade’s Population II stars. Thus, since the relative number 

‘ , Vy(x ; , ' 
of stars per unit volume of a given mass a =e-™(Y-)), the light emitted per unit 
Vy 
volume becomes 


a * 2 
kL 2 v(x) du= | wy (o)e™UDdu; v,(o)=p"%e-™. 


0 


Thus the light is proportional to 
= 00 
I. pes un-te ney du, i.e. to y iat 2) 


The solution of the equations gives the light and density distributions in space and, 
for comparison with observations, the projections of these against the celestial 


* D. M. Myers and W. R. Blunden, Fourn. I.E. Aust., 24, 10-11, 195, 1952. 
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sphere must be found. The projection of density is of the form | y "d(x? —a?)!2, 


where a is the projected central distance. The light projection is similarly 
«2 
| es +5/2) d (a* a a*)t?2, 
Ya 
It was found, by a comparison with the numerical solution for case n = 3, that 
the asymptotic solution differed from it by less than 0-2 per cent atx=20. For 
a>20, therefore, the projections were computed from 


( n—I ) 
2(n+1 
P, = pulate) q-@-Din+D (77/28 ey for the density 


n 
(55) 
I+” 
and n+4 
al 2(n+ 1) : 
P, = pian +5yi(an +2) q-(n+4)(n+1) (x/2)"? or for the light. 


r(= 5) 
2n+2 


Anticipating a little, another integral which it was found necessary to compute 





was . (x? — a?) dy-" with x»>a,and a>z2o0. This range of x enables us to use 
the asymptotic form of y, giving 

[et atprdy-n—yorm (apt —arye— |" y-nd (tat 
the latter negra being given by . 











- 2 
| . y-” d(x? —a2)'2 = pula 1m FP wat My an sayin 
ins I—n 3n+1 
4 
na (2n + 1) Xe +3)/(n+1) ‘ (3) 
2(n+1)(5”+ 3) 


A fuller discussion of the computational procedures is given in the following 
section. id 
It can be seen that the total light | P,,xdx and a fortiori the total number of 


a) 

stars | P, x dx in each of these models is infinite ; thus in the comparison of these 
0 

models with observational data, a method must be found to ascribe to them a 


finite radius. (Even so, as the projections themselves are finite, comparisons with 
observation can be made in the inner parts of the models as if there were a finite 
radius.) Truncation of the Maxwellian velocity distribution is fully equivalent 
to this both in the physical and analytical aspects of the problem. 

2. Solutions of the equations.—The solution was carried out in three parts, 
namely (a) determination of y in terms of x for various values of ; (6) evaluation 
of the density projections for various values of a and ; (c) evaluation of the light 
projections for various values of a and n. 

(a) There is no known analytical solution to equation (1) and a numerical 
solution for an adequate range of x and for all required values of m would be 
extremely tedious. It was therefore decided to obtain a set of solutions on the 
C.S.1.R.O. differential analyser. 

It was appreciated that (i) the wide range of y-" makes the equation unsuitable 
for solution in its present form, and (ii) the singularity at x =o makes it necessary 
to use an approximate numerical solution, e.g. in series, to get away from the 
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singularity. The difficulty (i) is easily overcome by putting u=Inx, w=Iny, 
and the equation reduces to 
++ = exp {2u—(n+1)w}. (4) 
It is clear that the solutions of equation (4) have as asymptotes straight lines 
whose gradients are given by w=2/(n+1), and substitution of this, with w=o, 
in (4) gives the equation of the asymptote as 
2u—(n+1)w=2(n+ 3)/(m+1)?. 
The machine solution started at u=o, x=1 in each case, the difficulty (ii) 
being overcome by solving (1) in series for y giving 








# nx* n(8n+5)x® n(122n*—183n+70)x® 
y=It+—- - 
3.3 8.3.4.9" ¥.6.7.88 9.9! 
‘. n(605n* — 1947n? + 1981n — 630)x!° 





g.1r!-... 
a result similar to that of von Zeipel.* 

The solutions on the differential analyser (see Fig. 2) show that for a small 
value of x (of the order of x = 20), the curves approach their asymptotes to within 
several parts in a thousand. As the overall accuracy of the differential analyser 
is of this order, the asymptotes may therefore be taken as the solutions of the 
equation beyond x=20 without loss of accuracy. This is of considerable im- 
portance in evaluating the density projections. 


Integrators 


| 2 3 4 







integrand 


MOTOR Argument 
Integral 





u 
w 

w 

-Jwdw 
2u-(n+i)w 


exp f2u-(n¢i)w) 









Sexp (n+!)w} du 


Adding 
Unit 


@ Denotes negative sign 





Adding 
Unit 


Gearboxes 






U 
OUTPUT 


Fic. 1.—Differential analyser set-up for equation w=— w— { tw dw+ {exp {2zu—(n+1)w}du. 
(Scale factors omitted.) 

The set-up of the differential analyser for solving equation (1) in the modified 
form given in equation (4) is shown in Fig. 1, in which the notation standardized 
by Bush? is used. It was convenient to integrate each side of the equation, 
so that it was put on the machine in the form: 


w= —w— foodw+ [exp {2u—(n+1)w}du. 
Integrator 3 is used to generate the exponential, its argument 2u—(n+ 1) 


being obtained from the independent and dependent variables through an adding 
unit. Otherwise, the set-up is quite straightforward. 


* H. von Zeipel, Kungl. Svenska Vetenskapsakademiens Handlingar, Band 51, No. 5. 
+ V. Bush, Yourn. Franklin Inst., 212, 4, 447, 1931. 
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The output table was used to plot —nw(=log,)y~") against u(=log,)x) as 
this gave a convenient scale for each curve and provided suitable input data for 
evaluating the density projections. These results are shown, together with the 
corresponding asymptotes, in Fig. 2, and are tabulated also, in the form of y 
against x, in Table I. 

TABLE I 
Machine solution of negative Emden polytropes of orders 
2, 2°5, 3, 3'5, 4 foro <x <20 


x ye Ya-5 ¥s 3-5 v4 
° 1'o 1'o I‘ 1'0 1'0 
0°25 10103 10103 10103 10103 10103 
o'5 10407 10405 1°0402 1 ‘0400 1'0397 
0-75 10891 1°0880 1:0869 1:0858 1°0847 
I'o 1°1522 1°1493 1°1463 1°1433 1'1406 
I'l 1183 1°178 1°1731 I‘170 1°169 
1°2 1°215 1°208 I‘2011 1°198 1°198 
1°3 1°247 1°239 1°2301 1°226 1'224 
1°4 1°280 1°272 1°2600 1°254 1'258 
1°5 1°313 1°306 1°2905 1283 1:287 
1°6 1°348 1°340 1°3216 1°310 1°314 
‘7 1°385 1°375 1°3531 1°338 1°340 
1°8 1°422 1°409 1°3849 1°365 1-366 
19 1°459 1°444 1°4170 1°393 1°392 
2'0 1°497 1°480 1°4494 1°422 1°417 
2°1 1°537 I°515 1°4819 1°451 1°444 
22 1°577 I°551 1°5143 1480 1°470 
2°3 1°617 1°587 1°5467 I°510 1°496 
2°4 1°657 1 ‘622 1°5790 1°540 1°523 
2°5 1°697 1-658 16113 1°570 1°549 
26 1°737 1°694 1°6432 1°599 1°575 
2°7 1-778 1°730 1°6749 1626 1601 
2°8 1819 1°765 1°7065 1°654 1627 
2°9 1°861 1°800 1°7378 1680 1°652 
3°0 1°903 1°835 1°7689 1°706 1678 
3°1 1°945 1-869 1°7997 1°732 1*703 
3°2 1-988 1°904 1°8302 1°758 1-728 
3°3 2°030 1°939 18605 1°784 1°753 
3°4 2°072 1°974 18905 1°810 1°778 
3°5 2°113 2-008 1°9202 1°837 1802 
3°6 2°155 2°042 1°9496 1°863 1°825 
$7 2°197 2'075 1°9788 1-888 1°847 
3°8 2°239 2°106 2°0077 I°QI2 1-868 
3°9 2°281 2°135 20363 1°936 1°889 
40 2°322 2°163 2:0646 1959 1°909 
4°2 2°404 2°239 2°1204 2-006 1°950 
4°4 2°483 2°305 2-176 2°053 1°995 
4°6 2°565 2°367 2°229 2°097 2°035 
48 2°646 2°432 2282 2°142 2°074 
5'0 2-718 2°493 2°334 2°186 2°115 
52 2°795 2°559 2°385 2°231 2°154 
5°4 2°873 2°620 2°435 2°274 2°190 
5°6 2°957 2°679 2°484 2°318 2°222 
5°8 3°035 2°733 2°532 2°359 2°257 
6-0 37111 2°791 2°579 2°396 2°292 
6-2 3°184 2°856 2°625 2°430 2°325 
64 3°261 2°916 2°670 2°466 2°356 


66 3°336 2°967 2°715 2°504 2°388 
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TaBLe I (cont.) 
x ye Vers Ys Ya-5 va 
68 3°409 3°013 2°759 2°542 2°421 
70 3°478 3°065 2°803 2°580 2°453 
72 3°550 3°116 2846 2613 2°483 
7°4 3623 3°182 2°888 2°647 2°510 
7°6 3°694 3°241 2°930 2°679 2°538 
78 3°762 3°293 2°971 2°714 2°566 
8-0 3829 3°341 3°012 2°750 2°594 
8-2 3896 3°393 3°052 2°784 2°622 
8-4 3°964 3°444 3°091 2°816 2‘650 
8-6 4°032 3°491 3°130 2°845 2°677 
8-8 4°100 3°538 3°168 2°875 2°703 
9°0 4°166 3°584 3°206 2°904 2°729 
9°2 4°231 3°632 3°244 2°931 2°755 
9°4 4°296 3°679 3°281 2°958 2°779 
9°6 4°361 3°728 3°317 2987 2°803 
9°8 4°426 3°780 3°353 3'018 2°827 
10°0 4°490 3°830 3°389 3°048 2°852 
II 4°801 4°055 3°562 3°200 2972 
12 5°105 4°276 3°727 3°336 3083 
13 5°419 4°485 3°884 3°460 3°189 
14 5°706 4680 4°035 3°570 3°289 
15 5986 4°870 4°180 3686 3°383 
16 6-248 5064 4°319 3°796 3°473 
17 6-519 5°256 4°454 3°907 3°562 
18 6-769 5°449 4°584 4009 3°644 
19 7027 5627 4711 4°113 3°725 
20 7260 5°777 4°835 4°217 3806 
yy ae T T T T T T T T 
OF 
2.0+ 
3.0r 
40 4 1 1 1 4 1 1 1 l 1 L 
0 0.2 04 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 
Fic. 2.—Machine solution of equation (4); (log y~")=w plotted against u=(log x). 
(6) The density integrals, defined by 
m2) 
Poa} _ (x¥/a*—1)!d(y"), (5) 


were obtained for various values of n and a. For a>20 the asymptotic solution 
was evaluated from the equation (2) and a machine integration was not required. 

















626 A. W. Rodgers and D. M. Myers Vol. 114 


For a<20, the evaluation of the integrals is a matter of simple numerical computa- 
tion, but for convenience the differential analyser was used, the curve already 
obtained for logy" being fed into the machine from an input table. One 
integrator was used to evaluate y-" from the input data, whilst another input 
table was used to feed into the machine the value of (x?/a? — 1)". 

However, it was necessary to decide on an upper limit of the integration to 
provide sufficient accuracy. Assuming that machine integration will proceed 
at least as far as the apparent coalescence of the input curve and its asymptote, 
the contribution to the infinite integral for values of x beyond the limit of the 
machine solution can be readily determined. 


This contribution is given by | (x? — a*)*d(y-") which is given above in 


Ly 
equation (3), where x, is the limit to which the machine integration is taken. 
The density integrals are shown in Fig. 3(a). 





0.75 T T T T T T T T T T T T 
ost - 


\ 
0.25m\— - 


0.0 «| 
fn=2 


7 n=2.5 








a a Se a a a ee a a ee ee 
9 202 © 6o 86 +. 2 40 60 #80 200 22:0 20 260 





z—> 


Fic. 3(a).—Logarithm of projected density in each model plotted against z. 


(c) The integral of the light projection is similar to that of the density projection 
except that y” is replaced by y-”, where m=n+8. As the light integrals 
converge more rapidly than the density integrals it is possible to simplify the 
procedure by considering the problem in two parts. 

(i) 1<a<z2o0. It has been found from actual computing that integration 
up to x=250, the maximum value of x for which the curves of y-” were plotted, 
leaves an error in the integral of less than o-o1 per cent, so that the machine solution 
can be accepted as sufficiently accurate. 

(ii) 20<a< oo. The plotted curves of y-™ do not differ perceptibly from 
their asymptotes when x 5 20, so that integration of the asymptotes gives a suffi- 
ciently accurate result without recourse to the differential analyser. 

The required integral is that given above 


Py=—a{_ (s#/at—1)!*ay™). (6) 


rma 
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The central values of the projections of light and density at a=o could not 
be found using the machine set-up described above (see equations (5) and (6)). 
Therefore, after having read off log y-" and log x from the machine output table, 
the results were transformed to give y as a function of m. Thus the integrals P,|%° 
and P,|<° were found by numerical quadrature and the remaining contributions 
P,|20 and P;|z, were found from the asymptotic solutions for y. The curves 
showing the light projections plotted against z are shown in Fig. 3(d). 
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Fic. 3(6).—Logarithm of light projection plotted against z. 


3. Conclusion.—It can be shown in Figs. 3(a) and 3(d) that although for a 
given set of parameters, e.g. total mass, central density and radius, the curves 
are rather different for different n, by suitable adjustment of central intensity 
and radius scale all can be made to agree over a wide range of radius with 
observations such as de Vaucouleurs’ empirical relation.* Thus the quest for 
more information regarding these parameters is necessary before a discriminating 
fit can be obtained. 

The results demonstrate that a differential analyser can be used satisfactorily 
for this class of investigation with a reasonable degree of accuracy. 
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Research School of Physical Sciences, Department of Electrical Engineering, 
Australian National University, University of Sydney. 
Commonwealth Observatory, 
Mount Stromlo : 
1954 August 31. 


* G. de Vaucouleurs, M.N., 113, 2, 1952. 








THE STABILITY OF A SPHERICALLY SYMMETRIC CLUSTER 
OF STARS DESCRIBING CIRCULAR ORBITS 


C. Gilbert 
(Received 1954 September 27) 


Summary 


Some general properties are derived for a stationary system of many 
gravitating masses of a type described by Einstein some years ago. The 
system is assumed to give the average field and the average density of 
matter for a globular cluster of stars. 

It is shown that a particularly smooth join of the field within the cluster 
to the exterior field in free space can be made when the density of matter 
tends to zero at the boundary. Also the condition for the stability of the 
masses in their orbits is shown to place an upper limit to their velocities, 
which in most cases restricts the gravitational mass of the cluster to a value 
which is less than its rest mass. 





1. Einstein (x1) has shown how to construct a stationary system of many 
gravitating masses, each describing a circular orbit under the influence of the 
gravitational field of all the masses. Such a system has spherical symmetry 
and may bear a resemblance to a globular cluster of stars in the sense that it 
gives a good approximation to the average density of matter and the average 
field for the systems occurring in nature. A particular solution of Einstein’s 
equations has recently been considered by Sen and Roy (2) from this point 
of view. 

The object of this paper is to consider, in greater detail than has been done 
by Einstein, the structure of such a cluster of particles, with particular reference 
to the stability of the particles in their orbits. 

The structure of a cluster depends primarily on the density of matter, which 
in genera! is a function of the distance from the centre. When this density 
function is known the gravitational field and the speed of the particles in their 
orbits can be determined. The derivation of expressions for these quantities is 
carried out in polar coordinates, which differ from the isotropic coordinates used 
by Einstein. These coordinates have a particular advantage over the isotropic 
coordinates when the density of matter tends to zero at the boundary of the 
cluster. It is shown in this case that, by suitable choice of the density function, 
the field within the cluster can be joined to the Schwarzschild field outside so. 
that the field tensors and their derivatives up to any order can be made 
continuous. ; 

The condition that each particle shall have a stable orbit places an upper 
limit to the velocities of the particles in the cluster, and it is shown that in most 
cases this limits the gravitational mass of the cluster to a value which is less 
than its rest mass. 

2. The material—energy tensor and the field equations.—We consider a system 
consisting of particles of rest mass m, each describing a circular orbit about the 
centre O of the cluster. The gravitational field of the cluster has spherical 
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symmetry about O and each particle describes its orbit under the influence of 
this field, the effect of collisions between particles being neglected. In order to 
have a continuous distribution of matter we make the proper number density 
of particles my oo and the rest mass of each particle mo, in such a manner 
that the proper density of matter py=mmnz, is finite at all points of the cluster. 

Since we are considering a static system the gravitational field within the 
cluster is given by 

ds? = — e* dr? —r? d0? —r* sin? 6 dd? + e” dt?, (1) 
where A and v are functions of r only, and r <a, r=a being the boundary of the 
cluster. 

Within a volume element at a point P of the cluster, the particles move in 
directions perpendicular to OP with the same speed V, and on account of the 
spherical symmetry of the system the same number of particles have their 
velocity vectors in any direction perpendicular to OP. Writing g,; for the 
fundamental tensor of (1), and x‘ (¢=1, 2, 3, 4) for the coordinates (r, 0, ¢, t) 
respectively, the speed of a particle having velocity components dx*/dt (« = 1, 2, 3) 
will be 


v= |(-20 =). (2) 


Also the material—energy tensor at P will be 


, dx’ dx* 
Ti = M8ik Fe? (3) 


where the barred expression gives the mean value of the energy and momentum 
of a particle for all the particles in the volume element at P. From (1) and (2) we 
find for a particle having a velocity V 
e'(dt[ds)?=(1— "V2, (4) 
and consequently dt/ds has the same value for each particle in the volume element. 
The mean value in (3) can now be evaluated and gives the components 
T3= T3= — 4poV*(dt/ds)?, (5) 
T$= poe"(dt/ds)? = poo; (6) 
where po, is the density of matter referred to axes fixed at P. The other 
components of 7% are zero. 
The components (5) and (6) must now be substituted in the field equations 
80 75= —G}+4Gg5, (7) 


in which units have been chosen which give the Einstein constant the value 87. 
Making use of the well-known components of (7) for the metric (1), as given 
in (3), we find in a straightforward manner 


v’=(e*—1)r, (8) 
a | I 
Sapw=e*(~ — 5) +a (9) 
e’V2=(e4—1)/2. (10) 
Since the square of the velocity of light in the transverse direction is e’, 


equation (10) gives the square of the ratio of the velocity of a particle to the 
velocity of light in the same direction. It can be verified from the equations of 
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the geodesics given in Section 4 that this value of V gives the velocity of a particle 
in a circular orbit of radius r. Elimination of e” from equations (4) and (10) gives 


A 
(5) =. (11) 
s 3-e 
The information required to discuss the structure of a cluster is contained in 
the equations (8)-(11). 
3. The gravitational field of the cluster.—We make the following assumptions 
concerning the constitution of the cluster, 
(i) The density of matter pop is finite, continuous and not zero at all points 
within the cluster, but may be zero at the boundary. 
(ii) The particles describe circular orbits about the centre of symmetry, with 
the exception of particles at the centre, which are at rest. 
(iii) The velocity of any particle is less than the velocity of light in the same 


direction. 
We define 


cr 
u(r) =47| poor? dr (12) 
0 
for all values of r. 
From (i), when r<a, p(r) is a continuous function of r tending to zero as fast 
as r?atr=o. Whenr >a, u(r) has a constant value (a). 


The interior field.—From (11), (ii) and (iii), 
1 <e*<3. (13) 
Also, since e* is finite at r= 0, (g) gives, after multiplication by r? and integration 
between r=o and r=r, 


etn — (14) 


r 





From (14) and (g), e*=1 and A’=o at r=o, and e*>1 for r>o. In the 
interval o <r <a, e’ is therefore a bounded positive function of r which may 
either be a monotonic increasing function of 7, or have at least one maximum 
value in the interval. 

We set A+ v=0 at r=a, to make the field continuous with the Schwarzschild 
solution at the boundary. Then from (8) and (14) we find 


T 
1 | dr +1n(1 — 2u(a)/a). (15) 
From (8) v’ =o atr=oand v’>o forr>o. From the boundary condition above, 
v is finite and negative at r=a. Integration of both sides of (8) between r=o 
and r =a then gives a finite negative value for v at r=o. 

From the above it follows that e” has a minimum value <1 at r=o and 
increases monotonically to a value <1 at r=a. 

The interior field is given by (12), (14) and (15) when the density pop is a 
known function of r. Also the discussion above has shown that the field will 
be free from singularities. 

The exterior field and the boundary conditions.—In order to take the most 
general static gravitational field for the empty region r>a, we assume that the 
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radial coordinate R is an arbitrary function f(r) of the coordinate r. The 
Schwarzschild solution may then be written in the form 
ds? = — e§ dR* — R*d6? — R* sin? 6 dd? + e"d??, (16) 
where 
eS=e’=1—2y/R, (17) 
» being the gravitational mass of the cluster. 

We write g,; for the fundamental tensor of (16) when ds is expressed in terms 
of the coordinates x! and their differentials dx‘, and we assume for the boundary 
conditions at r=a that the g,;, for the interior and exterior fields, as well-as their 
first derivatives, are continuous. 

For the goo and g3, components these conditions give 


f(a) =a, f'(@)=1, (18) 
and because we have set A+v=o at r=a, the g,, and g,, components will be 
continuous if »=y(a). Their derivatives will also be continuous if 


d 
v= Fra, v= (19) 


at r=a, and from (8), (9), (17) and (18) these conditions will be satisfied if the 
density at r=a has the value 
Po=f" (a)e*/4zra. 

If the density vanishes at the boundary of the cluster we must have f(a) =o. 
This case can be obtained by taking R=r. It can then be verified that in the 
case when the density of matter tends to zero at the boundary of the cluster both 
the interior and the exterior fields are given by the same expressions (14) and (15), 
the value of ,.(r) for the exterior field having the constant value (a) in accordance 
with (12). 

By forming the first m derivatives of the equations (14) and (15) it is clear 
that in the case poy =0 at r=a, the g,; and their derivatives up to the mth order 
will be continuous provided p”(r), p’’(r)...u"(r) are all zero at r=a. This 
will be the case if dpgo/dr, d®po9/dr?. ..d"-*po9/dr"-1 are all zero at r=a. In 
particular the Riemann—Christoffel tensor, which depends only on g,; and their 
first and second derivatives, will be continuous at r=a, if dpyo/dr=o at r=a. 

The exterior field also will be free from singularities provided .<a/2, and 
this condition is satisfied because e’ satisfies the conditions (13) at r=a. 

4. The condition for the stability of the particles in the circular orbits.—To 
investigate the stability of the particles describing circular orbits, we use the 
first integrals of the equations for the geodesics. For orbits in the plane @=7/2 
these can be written (3) 


e# (Z) +r (Fy -¢ (a) +1=0, (20) 
(21) 


(22) 


where h and k are constants of integration. 
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Substituting from (21) and (22) in (20) we find 
dr 


Aa ® 2.—% h* 
e (3) =—1+h*e"’— = = F(r,h,k). (23) 


For a circular orbit of radius r, F(r,h,k) and F’(r,h,k) are both zero. Let the 
particle describing the orbit be subjected to a small perturbation which causes 
it to follow a geodesic path having constants H and K which differ slightly from 
h and k respectively. In general there will be a value r=), differing slightly 
from the value for the circular orbit, such that 

F'(b, H, K)= —v'e~’K? + 2H?/B =o. (24) 
F(b, H, K) will be a small positive quantity and for r=5+.x, where x is a small 
displacement, 
F(r, H, K) = F(b, H, K)+ F'(b, H, K)(x?/2). 

The particle will oscillate radially near the circle of radius b if F’(b, H, K)<o, 

that is, if 
(—v"+v")e-"K? — 6H?/b* <o. 
Substituting from (24) this gives the condition 
(—v" + v'?—3v’/b)e~’<o. 

A sufficient condition for the stability of the system is therefore 

v" —v'? + 3v'/r>0. (25) 
The case F’(r,h,k)=o would require further investigation, otherwise (25) 
would also be the necessary condition for stability. 

It is possible that equation (24) cannot be satisfied for a value r=65. This 
will be the case if 

e-"=A+B/r', 
where A and B are constants. In this case the equilibrium will be unstable 
because (23) then gives only one positive value of r for which dr/dt=o. This 
value of v makes the left-hand side of (25) identically zero and does not therefore 
affect the sufficient condition for stability. 

From (8) and (25) the sufficient condition for stability can be stated in the 
form 

r\'e* > —(e*—1)(3 —e’). (26) 
From (9) and (26) we obtain the condition 
(e*— 1)(2e*— 3) <8zppr7e™*. 


In a region of space where py, =0 we find, from (5) and (6), Ti=o. It follows 

that in free space the condition for stability in a circular orbit is that 
e* <3. (27) 

From (10) and (27) we deduce that the velocity of a particle describing a stable 
orbit in free space must be less than half the velocity of light. On account of 
the continuity of the field at r=a, the condition (27) must hold at the boundary 
of the cluster in order to ensure stability for displacements which make r>a. 

We can also take (27) to be a sufficient condition for stability within the 
cluster, for then the condition (26) will certainly be satisfied. Since A’=o when 
e* has a maximum value, we then obtain from (g) and (27) the condition 


Strpo9 < 477?. (28) 
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For a particle at the centre of the cluster h =o, and the stability of the particle 
can be shown, from (23), to be a consequence of e” having a minimum value. 
5. Calculation of the proper mass of the cluster.—The rest mass of the cluster is 


M= |pod¥, (29) 


where d= is the proper volume element at a point P for an observer moving with 
speed V’, and the integration is taken over the whole cluster. 
For the metric (1) 


i= ¢ el+"2 sin 0 dr dO dé, (30) 


and dt/ds can be found by eliminating V? between equations (10) and (11). 
Substituting from (6) and (30) in (29) we find 


M= - | | Poot? {he%(3 —e*)}? sin 6 dr dO dé 
0 0 


: 0 


r@ du(r 
= [HO eye -pypear, 
From the definition (12) du(r)/dr is always positive, and therefore, when the 
condition for stability (27) is satisfied, we find that 
M-p 
o< ——— <0°057. (31) 
M 

The gravitational mass of the cluster is then less than its proper mass by less 
than 6 per cent of the rest mass. On the other hand, when the only restrictions 
on e’ are those given by (13), can have values which are greater than M. 

6. Discussion of some special cases.—It will be convenient to distinguish two 
classes of clusters, called the A and B clusters, which differ in the behaviour 
of e’ within the cluster. 

For the A clusters we assume that e’ is a monotonic increasing function of r. 
Since A’ >o the condition for stability (26) is satisfied within the cluster provided 
e’<3. The condition for stability at the boundary, however, gives the condition 
e’ <3, according to (27). From (g) at r=a, pog>o and the proper density py 
cannot be zero. From (10) the velocities of the particles increase with increasing 
distance from the centre to a value which is less than half the velocity of light 
at r=a. Also, according to (31) the gravitational mass of the cluster is always 
less than the rest mass of the cluster. 

A short calculation has shown that the special cases considered by Einstein 
and by Sen and Roy are both examples of A clusters, and are therefore stable 
systems. The restriction placed on the velocities of the particles by the condition 
for stability, however, restricts these cases to the configurations for which » <M. 

For the B clusters we assume that e’ rises sharply to a maximum value <3 
near r=o and then decreases to the boundary, in such a manner that pop, >o 
for r<a, and popy=0 at r=a. For these clusters the proper density of matter is 
zero at the boundary, and the particles having the highest velocities are con- 
centrated near the centre of the cluster. Since the condition (27) is satisfied, 
these are stable clusters, and according to (31) the gravitational mass is always 
less than the rest mass of the cluster, 


43 
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From counts of stellar images on the photographic plates of globular clusters, 
it is inferred that the number density of stars in general decreases from the centre 
outwards at first rapidly and then falling gradually to zero near the boundary (4). 
It would therefore appear that the best approximation to the field for these clusters 
would be obtained from the B rather than from the A clusters. In addition the 
B clusters enable a smoother join of the two fields to be made at the boundary 
of the cluster. 


1 North House, 
Longbenton, 
Newcastle upon Tyne : 
1954 September 25. 
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NOTE ON SPECTRA OF K-TYPE DWARFS IN THE REGION 
OF H AND K 


M. W. Feast and A. D. Thackeray 
(Received 1954 November 18) 


Summary 
The line Si 3905 is considerably strengthened in K-type dwarfs compared 
with giants. 
Emission cores to H and K are recorded on dense exposures as apparently 
a normal feature in dwarfs of type K2 and later. 





The southern stars « Ind (4"-7, dK5) and « Cen B (1™-7, dK1) have no such 
bright counterparts among K-type dwarfs in the northern hemisphere. With 
many bright giants available, their spectra offer the best opportunity for studying 
absolute magnitude effects at extremes of surface gravity. This note presents 
two results of interest which have been discovered in the violet region near 
CattK on spectra obtained with the two-prism Cassegrain spectrograph 
attached to the 74-inch Radcliffe reflector. 

The line Si3905.—Comparison of the dwarf «Ind and giant « Tau (both 
Ks) in the region near K at a dispersion of 18 A/mm showed immediately one 
outstanding absolute magnitude effect; the line Si 3905 is considerably stronger 
in the dwarf (Plate 2a, 6). No other line in the region shows a comparable 
change, the effect being enhanced by the fact that most lines tend to be somewhat 
stronger in the giant. The line is well separated from the adjacent line Fe 3906. 
Wave-lengths measured in the two stars and in a supergiant are as follows: 


e Ind (dKs) x Tau (gK5) e Peg (cKo) Identification 


3905 °51 05°56 05°53 Sil 05°528 
3904°78 04°72 04°76 Til 04°785 
390388 03°95 03°84 Fe I 03:900 


Microphotometer tracings and visual examination suggest the possible 
presence of a blending line about 0-2 A to the violet of Si 3905 in the supergiant. 
We have not been able to identify this feature by consulting laboratory wave- 
lengths or Miss Davis’s measures in 8 Peg and «Sco (1). In «Ind and «Tau 
at least there is good reason to suppose that the measured line is predominantly 
due to silicon with negligible blending at the core. Photometric measures of 
the central intensity of Si3905 are as follows: 


e Ind x Tau e Peg 
0°233 0°432 0°356 


The tracings are unsuitable for measurement of equivalent widths owing to 
blending in the wings, particularly in the supergiant. Nevertheless it may be 
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estimated that the ratio of equivalent widths in dwarf and giant is ap proximately 
15:1. Further measures should be made when coudé dispersion becomes 
available. The strengthening of Si3g905 is not a peculiarity of «Ind. It has 
been confirmed in many other dwarfs by visual inspection. 

The most complete comparison of spectra of giant (Ko) and dwarf (Kr) 
spectra at high dispersion has been undertaken by Miss van Dijke (2). It 
was found in general that the lines below about 2-6 volts E.P. are strengthened 
in the giant, while lines of higher excitation are strengthened in the dwarf. 
The present work shows that Si3g905 (E.P. 1-9 volts) clearly violates this 
empirical rule. No lines due to silicon were included in Miss van Dijke’s 
comparison. Application of Miss van Dijke’s values for temperature and 
pressure in giant and dwarf shows that silicon should be predominantly neutral 
in both atmospheres; this conclusion should hold good for the cooler stars 
considered here. ‘The ionization potentials of Si (8-11 eV) and Fe (7-86 eV) 
are so similar that the marked difference of behaviour of Si 3905 and neighbour- 
ing Fe lines of low E.P. constitutes a departure from the predictions of ionization 
theory on the basis of constant abundances of these elements. 

A similar situation exists in the case of Mg 5183 which, as previously shown, 
is markedly stronger in K dwarfs (3). 

While this absolute magnitude effect shown by Si 3905 is thus of interest 
in its bearing on ionization theory, it is not likely to be of practical importance 
as a luminosity criterion since on low dispersion the line would become severely 
blended; and, further, it lies in a region of weak radiation. The strengthening 
of Mg 5183, 5172, 5169 in K dwarfs could however be a useful luminosity 
criterion on low dispersion. 

One possible explanation of the weakness of Si 3905 in giants is that the line 
may be partly filled in by emission. ‘The line is a prominent emission feature 
in Me variables and the other lines Fe 4202, 4308, well known as belonging to 
this category, have been shown by Spitzer (4) to have their absorption cores 
partially filled in by emission in the supergiant «Ori. This explanation can 
be tested only with higher resolving power than is at present available at the 
Radcliffe Observatory. 


TABLE I 
Star Sp. M, Ca II emission 

pEriA dK2 6-9 doubtful 

p Eri B dK2 70 certain 
*e Eri dK2 6-2 certain 
*§ Eri dKo 3°8 no 
*40 Eri A dK1 6-0 doubtful 

y Lep B dkK3 6-7 certain 

A Vel cKs certain 

Proxima Cen dM 15°7 certain 

a Cen B dK1 6:1 certain 

HR 7703 A dK2 6°5 doubtful 
*61 Cyg A dK3 74 certain 
*61 Cyg B dKs 8-4 certain 
*e Peg cKo certain 

e Ind dKs 70 certain 
*y Tau gKs certain 


The spectral types and absolute magnitudes of the dwarfs are from Kuiper (Ap. 7., 
95, 201, 1942). ‘The Yerkes type for ¢€ Peg is K2Ib. The stars marked by an asterisk 
appear in (5). 
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H and K emission.—It is well known that many late-type stars, dwarfs and 
giants, show emission cores to the Catt H and K lines. In particular, Joy and 
Wilson (5) have listed 445 stars showing such emission lines. "The phenomenon 
becomes increasingly common among the latest types, especially dwarfs. Joy’s 
plates for radial velocity showed that the phenomenon became general at type 
dM4°5 and later (6). 

The Radcliffe plates taken to study the behaviour of Si 3905 are sufficiently 
exposed to show that H and K emission is a normal feature of dwarfs of type K2 
or later. ‘Table I presents our findings concerning Cali emission in 15 stars. 

It is noteworthy that emission was not certainly recorded in p Eri A on the 
same plate which clearly recorded it in p Eri B, although the two components 
have practically the same brightness (Plate 2e, f). Plates of « Ind are available 
covering two seasons but have not as yet shown any sure signs of variation in 
the intensity of the emission. 

It should be emphasized that the Ca11 emission is only recorded by dint of 
overexposing for the normal region (Hy to Hd). Given adequate exposure 
and resolving power the emiss‘on cores could probably be detected in normal 
stars of still earlier type. By analogy with the Sun the emission probably 
arises at a chromospheric level. 


Radcliffe Observatory, 
Pretoria : 
1954 November 8. 
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ELECTROMAGNETIC FIELD EQUATIONS 
FOR A MOVING MEDIUM WITH HALL CONDUCTIVITY 


F. H. Piddington 


(Received 1954 August 17) 


Summary 
In an isotropically conducting medium the electric field and current 
density are associated by Ohm’s law. This, together with Maxwell’s 
equations, leads to - 


oH 7H 
V?H=4704 — —curl (Ux H)>4+-c"* => > 
{ ot } ot* 
where H, v and o are the magnetic field, velocity of the medium and con- 
ductivity. ‘This equation describes any electromagnetic field in a moving 
(or fixed) conductor. It includes the “‘ telegraph equation ”’ of radio waves 
and also describes any hydromagnetic wave. 

Ionized gas in a sufficiently strong magnetic field is anisotropically 
conducting and o is replaced by a tensor. An equation analogous to the 
above is derived for these conditions, allowing the study of electromagnetic 
disturbances in a medium with Hall conductivity. The equation is used 
on some simple problems of astrophysical interest, including hydromagnetic 
waves. 





1. Introduction.—The dynamics of an ionized medium may be studied by 
considering the microscopic motions of individual particles and summing their 
effects. When a strong magnetic field is present, however, the resulting equations 
are so complicated as to be insoluble except in some very simple cases. 

Numerous more complicated problems are met, particularly in the field of 
astrophysics, and the most promising approach to these appears to be in the use 
of “‘transport equations’’. ‘These statistically combine individual motions and 
describe gas movement on a macroscopic scale. ‘The transport equation used 
here is that relating electric current density j to the impressed electric field E. 
In its simplest form, when the conductivity is a scalar quantity o, when the 
medium is at rest and when E does not vary rapidly, it is simply Ohm’s law 

j=cE. (1) 
When the medium moves with velocity v in the presence of a magnetic field H, 
it takes the form (in electromagnetic units) 


j=o(E+v~ H). (2) 


~ 


‘The current so determined distorts the magnetic field according to the 
equation 


54 20k 
curl H= 47j +c a ( 


Combining (2) and (3) with Maxwell’s other equations 


w 
~— 


curlE=~ 5" and divH=o, (4) 
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an equation is found describing the rates of change of the electro-magnetic 
(including radiation) field 
CH 


V7H = gral —curl(v x #)} +c? ap ° (5) 


This well-known equation, or one of its simpler forms, is the key to a wide 
range of problems which fall into two categories. (a)—When v is zero or has 
prescribed values which are independent of the electromagnetic field, if the 
last term is important it becomes the telegraph equation describing radio waves 
in conductors which are fixed or moving in a known manner; if the last term 
is negligible it describes the diffusion of a magnetic field in a rigid conducting 
medium. (b)—When v depends appreciably on the electromagnetic field, the 
current j, as well as influencing the field, also causes a mechanical reaction j x H 
on the medium which appears in the hydrodynamic equations. ‘These, together 
with equation (5), then describe any magneto-hydrodynamic (or hydromagnetic) 
phenomenon including the type of wave motion first investigated by Alfvén. 

Unfortunately this relatively simple situation does not obtain when the 
magnetic field becomes so powerful that the radius of gyration of the conducting 
electrons or ions is of the same order as (or smaller than) their mean collision 
distance. The medium becomes anisotropic, the conductivity being a tensor, 
and equation (1) is replaced by 


. E 
j= )E,+0,E,+0,H x a (6) 


H? 


where E, and E, are the components of E parallel and perpendicular to H. 
The conductivities oy, o, and o, give rise to current components along H, 
parallel to E, and perpendicular to both H and E, (Hall current). 

As for the simple medium, equation (6) is one of a chain of equations 
linking j, H, v and E and, in principle at least, allowing all the variables to be 
simultaneously determined. One approach to this difficult problem is to neglect 
the space derivative of the magnetic field caused by the current j. ‘Then the 
induced electric field given by equation (4) and the consequent change in jf 
given by equation (6) are both neglected, thus greatly simplifying the problem. 
This method, valid of course only under certain conditions, has been used 
elsewhere (1) to find v, j and the internal electric field in a medium subjected 
to external magnetic, electric and mechanical fields of force. ‘The other approach 
to the problem which does include consideration of induction effects is to derive 
an equation analogous to (5) and capable of being used like (5). 

In the present paper equation (6) is combined with Maxwell’s equations 
to obtain the electromagnetic field equation for a moving, anisotropically 
conducting medium. The equation resembles (5) and in many problems reduces 
to an equally simple and almost identical form. 

The results are applied to some astrophysical problems including hydro- 
magnetic waves* for which the rate of dissipation of energy may be determined. 
The results may also be usefully applied to some non-oscillatory problems. 

* In private discussion, Dr D. F. Martyn has reported an investigation of wave motions in 
parts of the ionosphere where Hall conductivity is important. His results seem in general agreement 


with those given below, referring to the particular geometrical configuration with which he is con- 
cerned, 
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These include further extension of the magnetic storm theory of Chapman and 
Ferraro, the movement of sunspot magnetic fields in the chromosphere and 
corona and perhaps the movement of ionized interstellar gas. 

Although the examples all refer to gases, the equations developed are equally 
applicable to liquids or solids and might, for example, have useful application 
to semi-conductors in strong magnetic fields. 

2. The electromagnetic field equations of a non-isotropic medium.—Consider a 
medium permeated by a steady, homogeneous magnetic field of external origin 
and strength, Hy, such that the Hall conductivity o, is appreciable compared 
with o,. The variable electromagnetic field also present has a magnetic field 
component much weaker than H, so that the total magnetic field H= HM. It 
is also assumed that the conductivity components are all constant within the 
region concerned and that changes in the electromagnetic field are slow enough 
to allow the use of the transport equation. These restrictions are discussed in 
Section 3. 

Choose a right-handed system of axes with the z axis lying along H, and 
use the subscripts x, y and z to define the components of the field along those 
three axes. It is desirable here to introduce a value of conductivity defined by 
the identity 0; = (¢,+,7/0,); this is discussed in Section 4. 

The current components are given by 


J2=%F,—o2F, 
Jy =F y+ oF 
J2=%F,, 
where F=E+v0x i. (10) 


Let us, for the time being, neglect the displacement current term in equation (3) 
so that 
curl H= 47j, (11) 
from which 
\V?H = — 47curlj. (12) 
The x components of this equation give 


ay O. Ay 
Viera ae & ty i) 


= = 40 (0528 =o, 0,222) - 
ee 70 Oy 9: "2 az J° 3 


Differentiate equation (7) with regard to z, multiply by 470,,c, and subtract 
from (13), 
2 Of OF, OF OF, 

v2 = >) el Fs os ae 

1H — 40 a, Oz gro, (5 ds ) 47(%— 03) dy © 
Using equations (4), (9) and (11) this reduces to 

2 a, 0 . ey? a3\ 0 

V?7H,- i Te gem, FY curl, (vx H)}—( — “2) a (14) 


% 


A similar equation in H, may be derived starting from the y components of 
equation (12). 
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Now consider the z components of (12) 


shitciner teeta sna aaah 


te  *% +o, y (15) 


Differentiate equation (7) with regard to x, (8) with regard to y and add. For 
the relatively slow changes in the electromagnetic field (see Section 3) the 
medium remains everywhere neutral, so that divj=o. Hence we have 


oF, rs oF, rs oF, _ oF, di oP = 
“1\ Ox oy “s oy ox az” 
which, when multiplied by 470,/0, and combined with (15) gives 


V?H,- Fe curl, H= gross —curl,(v x H)| . (16) 
o, Oz ot J 

Equations (14), (16) and the corresponding one in H,,, together give the full 
vector equation of the rates of change of the electromagnetic field. ‘Their 
generality is subject only to the restrictions mentioned above. However, the 
additional terms to those in equation (5) (space derivatives of the components 
of curl H), may add greatly to the difficulty of solving the equation and it is 
desirable to examine further restrictions which provide simplification. 

In the first place, when dealing with a gas for which only electron conductivity 
is important, 0, oy. Such a gas, although it may include substantial pro- 
portions of other ions, has the approximate characteristics of the binary gas 
discussed in Section 4. It is a gas frequently met in astrophysical problems. 
The equations are simplified for such a gas by omission of the term in (1 —a3/a,). 

For other types of conductor the same terms are omitted, provided the space 
derivatives of curl, H may be neglected. But curl, H indicates a twisting of the 
magnetic lines of force about the direction of the main part of the magnetic 
field H). Thus the term in question indicates a variation in the xy plane of such 
a twist. There is no doubt that such effects occur and they may be important, 
in which case the full equation must be used. In other problems, or for a binary 
type gas, the field equations reduce to the form 


\7?H — 02/0, < curl H= 470, oS —curl(v x #} ; (17) 


Finally, for problems where d/dz=0, that is for disturbances propagated 
across the magnetic field, the second term in equation (17) may be omitted. 
When this is done the displacement current term may, if desired, be reintroduced 
as in equation (5) and we have 

VIM = 404 {SF —curl(vx Hb +e #2 Z (18) 
at ot? 
This is of the same form as equation (5), merely having co; substituted for a, so 
that problems of this nature have similar solutions to those for an isotropic medium 
and may be solved by identical procedures. All of the examples given below, 
admittedly geometrically simple, are solved by the use of equations (17) and (18). 

3. Discussion of the equations.—A restriction on the use of any of the equations 
derived above is that the external magnetic field must be steady, homogeneous 
and large compared with the variable component. This restriction is necessary 
because the space and time derivatives of the conductivity components have been 
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neglected. ‘hese components depend on the magnitude and direction of the 
total magnetic field. It is necessary, therefore, that the derivatives of H/H, 
should be small compared with the corresponding derivatives of the electro- 
magnetic field components. 

Since o, and o, also depend on the density and temperature of the various 
components of the gas (see Section 4), these must also be stable within the same 
limits. 

These conditions are usually met for oscillatory disturbances, such as hydro- 
magnetic waves of small amplitude or radio waves, because the wave-length is 
much smaller than the extent of the medium and of the external magnetic field. 
There is another class of problem, however, in which the gradients of the 
‘‘external’”’ magnetic field (and perhaps also the gas parameters) are of the same 
order as those for the electromagnetic field. Such a problem is the large-scale 
movement of sunspot magnetic fields and their associated gas in the corona. 
The magnetic field is no longer external but part of the variable electromagnetic 
field. ‘The equations of Section 2 are not applicable to these problems. For- 
tunately, however, most such problems met in astrophysics are concerned with 
a gas in which only electron conductivity is important (see Section 4). For such 
gases o, depends only on gas temperature and if this can be assumed uniform 
then many problems may be solved as shown in Section 7. 

The values of o, and o, are derived on the assumption that the imposed 
electric and magnetic fields do not vary appreciably in a time necessary to 
provide a good statistical average of the ion motions. When equation (1) is 
relevant the time concerned is large compared with the ion collision periods. 
When the ions and electrons complete one or more rotations in the magnetic 
field between each collision it is the gyro-period which determines the averaging 
time necessary. ‘The latter must be long compared with the gyro-period of each 
ion which contributes appreciably to the conductivity. In the case of a lightly 
ionized gas (see following section) both electrons and heavy ions may contribute 
to the conductivity, so that the time must be long compared with the gyro-period 
of the latter. In a binary gas only electrons contribute appreciably to the 
conductivity, so that their gyro-period is the deciding factor. 

In very tenuous gases another effect may be important in limiting permissible 
rates of change of the variables. In deriving equation (16) it was assumed that 
divj=o, which is not strictly true. Disturbances (let us call them waves) 
propagated across or obliquely to the magnetic field are associated with plasma 
waves. While the latter may be important and their electric fields substantial, 
the current required to maintain them is usually insignificant, in fact comparable 
with displacement current. In the present paper this is assumed so, but criteria 
for determining the effect of the plasma wave current are given in a paper in 
preparation. 

4. The conductivity.—The most complete and detailed account of the deriva- 
tion of conductivity (and other) transport equations from the Boltzmann equations 
of gas mixtures is given by Chapman and Cowling (2) and Cowling (3). With 
their values of o, and o, the value of o, is found from the identity given in 
Section 2. Its physical significance was first pointed out by Cowling (4)*, who 

* While Cowling’s proof was not restricted, his paper dealt with a heavily ionized gas. The 


significance of o; in the case of a lightly ionized gas (when o3>409) may first have been realized 


by Baker and Martyn (5). 
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showed that when Hall current is prevented from flowing, then the conductivity 
parallel to E, becomes o; instead of 0}. 

For convenience the formulae for og, o,, o, and o, for a binary gas and a 
lightly ionized gas have been summarized in reference (1). In general they all 
depend on the number density of the gas components, on the gas temperature 
and on Hp), but in the case of a binary gas o,=0, and depends only on gas 
temperature. This important result is used in Section 7. 

In Table I a list of conductivity components is given for some regions of 
astrophysical interest specified by the electron density n, and gas temperature 7’. 
In interstellar space and most of the solar atmosphere the binary gas formulae 


TABLE I 


Conductivities in some regions of astrophysical interest 





Region 
Gas 
parameters 





Interstellar 
H I region 


i 
Interstellar | 


H II region 


Solar 
corona 


Solar 
chromosphere 


Ionospheric 


region E 


lonospheric 
region F2 


(100 km) 








10° 10° 
300 10° 


03 03 


10-* 10 10° 
100 10! 10° 
s07* 10-3 30 
KX 107-12 x 107° 6 x 10-8 : x 107} >1071 
K yo" x r3o~** 6x 107* z0~** 4x z0-*9 
<1o717 < 19716 610714 -13 7 10719 16 
707" K 207? 6 x 1078 : 10 axa 























hold closely enough and in the magnetic fields given we have o,( =c,)) >o,> 0}. 
‘This means that when equation (18) is applicable the results are the same as 
for an isotropically conducting medium of conductivity o,; that is, the same 
medium in the absence of a magnetic field. Also the “effective” conductivity o, 
is much larger than either o, or og. As the value of H, is reduced so that the 
electron gyro-frequency approaches the collision frequency, o, increases con- 
tinuously and o, increases to a maximum and then decreases. When gyro-angular 
frequency and collision frequency are equal, og=0, and o,=20,; also oy has 
its maximum value. In yet weaker magnetic fields o,/0, becomes small and 
the main contribution to a; =o, is made by o,. 

The position is different for a lightly ionized gas, as shown in the last two 
columns of Table I, calculated by Baker and Martyn (5). When the gyro- 
frequency is large compared with the collision frequency (F2 region), 0, <0, <ay, 
so that o,~o,. This marked difference in the behaviour of tenuous binary and 
lightly ionized gases is discussed in reference (1). For the latter, equation (6) 
loses its last term and equations (14), (16) and (17) are simplified by losing the 
terms in o9/0. 

In the following sections the equations of Section 2 and the data of ‘Table I 
are used to discuss some problems of electrodynamic disturbances. Such 
disturbances occur whenever the electromagnetic field and forces acting on 
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the gas are not in equilibrium, either static or dynamic. ‘They may be roughly 
divided into small-scale and large-scale disturbances, referring to the extent 
and not the violence of the disturbance. A small-scale disturbance is one in 
which the main variations in the electromagnetic field occur in regions small 
compared with the total extent of the gas and of the steady magnetic field. 
Hydromagnetic waves usually fall into this category. A large-scale disturbance 
extends throughout a substantial portion of the magnetic field which now cannot 
truly be described as either “external” or “steady”. The additional difficulties 
met in such problems are discussed in Section 7. 

5. Disturbances propagated along the magnetic field.—Consider a transverse, 
plane wave (either a single disturbance or train of waves) propagated along the 
magnetic field H,; that is, along the z axis. ‘Then, since 0/0x=0/dy=0, 
equation (17) may be used. Neglecting second-order terms the x and y 
components become 














,@H, OH, 4, 00, _ so, PH, ca 
° Oat ~ Ot 7° Gs ~ o, Oat” 9 
@H, OH, ,, dv, , sop PH, 

TA ee ee, Oe: (20) 


When the conducting medium does not move (v,=v,=0), because of adequate 
rigidity or inertia, the equations reduce to a modified “telegraph equation”’. 
For oscillatory disturbances they indicate two waves, circularly polarized in 
opposite senses. ‘The dispersion equation for solutions of the form exp i(pt — kx) 
is 


al 4703 


~ p( +o9/0,+i) 


| > 
tre tw 


This equation should show agreement with the corresponding magneto- 
ionic equation. Substituting values for the conductivities of a gas having only 
appreciable electron conductivity (the ions are assumed immobile), 


k2¢2 — 47n,e7c? 
p  — mp( Fw +i/r) * 


where é, m, w and 7 are the electron charge (e.m.u.), mass, gyro-magnetic frequency 
and collision period and the other terms have been defined. Agreement is shown 
with the magneto-ionic theory except for the omission of two terms on the r.h.s. 
A term unity is missing due to the neglect of the displacement current and a 
term p inside the brackets is missing because p has been assumed negligibly 
small compared with w. The effect of displacement current is not that described 
in Section 3 above, which occurs only when the gas is compressed (oblique or 
transverse propagation) ; it is usually negligible even for frequencies approaching 
the ion gyro-frequency. 

When the medium is moved by the electromagnetic forces, hydromagnetic 
waves of the kind first investigated by Alfvén (6) may occur. Neglecting 
compression the equations of motion of the medium are 


dv, V?0H, dv, V?0H, 


and —! 





ot =H, os ot H, dz’ 
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where V?=H,?/47p and p is the mass density. Solving with equations (19) 
and (20) two hydromagnetic waves are found, circularly polarized in opposite 
directions. The dispersion equation is 


o (ves sp | + i 


p? ~ 470,90; 47705 


In the case of a binary gas (to which the solar corona approximates closely 
enough) the ratio of the first two terms within the brackets is (2/p, where Q is the 
heavy-ion gyromagnetic frequency. Hence the ratio is very large and the two 
waves have nearly equal velocities V +o,p/870,0,V or V + H,/4An,e, where A 
is the wave-length. The polarization and velocity of these waves are the same 
as those given previously by Cowling (reference (9), p. 582). 

In the case of a lightly ionized gas the velocities may differ widely and when 
V? <oyp/470,0, one wave will not be propagated. It should be noted that such 
waves require that each neutral molecule should experience many collisions in 
one wave period. If many wave periods occur within the collision period of a 
neutral molecule then the hydromagnetic waves are propagated in the electron-ion 
plasma and not in the gas as a whole. 

In all cases the absorption coefficient is p?/870,V, thus depending only on o, 
among the conductivities. 

This result is of considerable significance in astrophysics. It was previously 
often assumed that o, was the conductivity mainly effective in determining the 
losses sustained by a hydromagnetic wave.* It is evident from Table I that o, 
is often very much larger (by factors up to 1o!*) than o,. Even had o, been 
taken as the relevant conductivity, the error would still be very large. To 
illustrate this point we refer to a theory of Alfvén (8) concerning the propagation 
of hydromagnetic waves from the solar photosphere to the corona. According 
to Alfvén such waves would be absorbed in the corona, their energy then being 
available to maintain the observed high coronal temperature. However, in 
calculating the absorption coefficient, Alfvén neglects Hall conductivity and 
assumes o, to be the determining conductivity. For the conditions he assumes 
(similar to those in Table I) o,>6x 10° e.m.u. while o, > 6x 10° e.m.u. 
Hence the rate of absorption is overestimated by a factor of ~10'*, the true 
rate being negligible. ‘This conclusion was anticipated by Cowling (9g) who 
pointed out that the Joule loss, which depends on collisions, could not be 
increased by the presence of a magnetic field. 

6. Disturbances propagated across the magnetic field—Consider a plane 
electromagnetic-hydrodynamic disturbance (including the case of a train of 
waves) to be propagated across the magnetic field H), say along the x axis, so that 
d/dy=0/dz=0. Such a disturbance is satisfactorily described by equation (16), 
which is simplified because the term in 0/dz curl, H is zero. Also 


curl, (vu x H)= — H, (dv,/0x), 
so that the equation becomes 


0H, OH. dv, 
oa = 05x" (21) 





* Schluter (7), studying binary gas by a combination of microscopic and macroscopic methods, 
concluded that og was the relevant conductivity, 
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This equation describes any electromagnetic disturbance in an anisotropically 
conducting medium propagated across the magnetic field and subject to the 
restrictions of Section 3. It includes the telegraph equation and the equation 
of hydromagnetic waves. Again o3, rather than o,, is the conductivity which 
determines the rate of dissipation of energy. 

In order to study the disturbance represented by equation (21) a second 
relationship between v, and H, is required and this is provided by the hydro- 
dynamical equations of momentum and continuity and the gas equation. These 
are 

dv, 00 H, 0H, 


Ue a 
Poe TP ea t ae = 47 Ox °* 





0 Op 
ox (pv,)+ 57 = 


0=f(p, T), 
where @ is the gas pressure, being a function of density and temperature. In 
general the equations are intractable, but for vanishingly small waves 0/dt > v( 0/0x) 
and f=const.p”, where y is the ratio of specific heats. The three hydrodynamic 
equations then combine to give 
OV, 77_0*V, V2 0°(H,*) a 
on dx? 2H, oxot (a8) 
where U?=06/dp is the square of the velocity of sound waves in the gas (in the 
absence of a magnetic field). Equations (21) and (22) may be solved simul- 


taneously by assuming a solution of the form expi(pt—kx). The dispersion 
equation is 














ws —fe{ pays POTN mo. (23) 
| sU? sU? 

This equation, not previously derived, allows the study of small transverse 
hydromagnetic waves in an anisotropically conducting medium. It cannot be 
compared directly with the corresponding magneto-ionic equation unless it is 
assumed that p<470,c*. This is due to the effect of displacement current dis- 
cussed in Section 3. It may be compared with an equation due to Anderson (10), 
derived for an isotropically conducting gas. The only difference is in the value 
of s which, for a scalar conductivity a, is (470)"}. 

For permissible waves in the gases of Table I (wave-length > 1 cm) we have 


p(U?+V?). pp’ 


> ca 


The roots of equation (23) are given approximately by 
p? _ ip(U*+ V*) p° 
i, a en nee 
_ U? sU? _ (U?+ V?)+ isp” 


The first pair of waves are rapidly damped, the second pair give 


k= +P I - 
~ (U?+ Vy 2(U?+ V?)J° 
This comprises a pair of waves travelling in opposite directions with velocity 
(U?+V*)"? and absorption coefficient p*s/2(U?+V?)**. Because the damping 
has here been assumed small, the expression for the velocity is the same as that 
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previously determined (11, 12) for a perfectly conducting gas. If the damping 
is large then the velocity will depart from this value. 

It is of interest to calculate the rates of attenuation of hydromagnetic waves 
propagated across the magnetic field for some gases of ‘Table I. In each case 
the wave frequency is chosen to be as large as possible but still an order of 
magnitude less than the ion gyro-frequency. In the chtomosphere the latter 
(for protons) is ~ 10° rad. sec"!, so that a wave of frequency 10° rad. sec ! is 
considered. Its velocity, mainly due to the resilience of the magnetic field 
(V>U), is 4x 108 cmsec™! and it diminishes by a factor exp(—1) in a distance 
of about 5x10'km or 2x 10° wave-lengths. Observe that the wave may be 
propagated only because of the high value of o, and would be damped within 
a wave-length if it depended on o,. Waves of lower frequency will suffer less 
damping. Waves propagated in the corona will not be damped appreciably in 
distances large compared with the linear dimensions of the corona. 

In the region listed in Table I as interstellar H II, the proton angular 
gyro-frequency is about 10 rad. sec!, so that hydromagnetic waves of the type 
described by equation (23) might have frequencies as high as, say, 1 rad. sec"! 
(wave-length ~ 4000 km). Such waves would travel about 0-01 parsec before 
being attenuated by a factor exp(—1). 

Hydromagnetic waves of length short compared with the thickness of 
ionospheric E or F layers are damped too heavily to allow propagation at all. 
Longer waves introduce boundary problems with which we have not concerned 
ourselves here. 

7. Large-scale disturbances.—A typical large-scale disturbance (according to 
our definition) occurs when a sunspot magnetic field and its associated gas enters 
the solar atmosphere. Here the gradients in the electromagnetic field parameters 
have the same order of magnitude as those of the total magnetic field and of the 
gas parameters. Hence one of the restrictions listed in Section 3 is not satisfied 
and the equations are not valid. Another difficulty usually met in these problems 
is the complicated geometry of the system of magnetic fields and gas distribution. 
These difficulties are usually insuperable, but an insight intothe probable outcome 
of the disturbances may often be made by assuming a simplified geometry. This 
method was used with considerable success by Chapman and Ferraro (13) in 
their investigation of a complicated large-scale disturbance: the initiation of a 
geomagnetic storm by a jet of solar gas. 

Let us first consider the simplest possible geometrical configuration of a 
disturbance, in which the electromagnetic and gas parameters vary in only one 
direction. Let this be perpendicular to the magnetic field, say along the «x axis. 
It is sufficient to consider the z component of equation (12): \/?H,= —47curl,j. 
Since 0/0y=0/0z=0 this reduces to (writing H for H,) 

. ; 
oa - — gn 2, (24) 


Now curl, H is zero and hence j, is zero and equation (7) gives F,,=(a,/0,)F,. 
Substituting in equation (8) we have 


Jy =F 


y’ 
which, with equation (24), gives 


07H he cio a 
i i Ox 


Ox 


(o,F',). 
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In the case of a lightly ionized gas the space derivatives of o, may be important, 
in which case a solution may be difficult or impossible to obtain. However, in 
most problems of astrophysical interest, the gas may be considered a binary gas 
and o, depends only on its temperature. Assuming isothermal conditions, 
equation (25) reduces to the well-known form 


oH OH Oa 
Fu = 4705 Es + =, (ct) } (26) 


The result indicates that transverse disturbances in an isothermal binary gas 
show induction effects depending only on the gas temperature. The value of 
the total magnetic field at any time and place, although it determines o, and og, 
does not affect the rate of diffusion of the electromagnetic field. It does, however, 
cause a physical effect not shown in equation (26) but discussed elsewhere 
(1, 14). The gas does not move across the magnetic field in the direction of the 
mechanical force causing its motion, but at an angle tan-'a,/o, to this force. 
Thus, while the motion of the gas is determined by the magnetic field strength, 
the ‘“‘motion”’ of the electromagnetic field is not affected. 

The next simplest type of disturbance after a plane wave is a cylindrical 
wave, the magnetic field being parallel to the axis of symmetry. This will not 
be considered separately but as a special case of a more complex system in which 
the magnetic field varies in direction as well as intensity. The restriction imposed 
is that the magnetic field is symmetrical about an axis (say the z axis) and the lines 
of force all liein planes through Oz. Such fields include those of dipoles, uniformly 
magnetized spheres and, with sufficient accuracy for illustrative purposes, sunspot 
and stellar fields. ‘The mechanical forces associated with such fields will also be 
symmetrical and lie in planes through Oz. The electric currents will flow in 
circles with centres on Oz and lying in planes perpendicular to Oz, so that the 
force j x H is everywhere in the meridian plane and perpendicular to H. There 
can be no accumulation of electric charge and hence no boundary effects caused 
by such accumulation. 

It is evident that the above electromagnetic field violates the first restriction 
applying to the equations of Section 2. The magnetic field varies not only in 
strength but in direction from place to place. Fresh equations must be derived 
which satisfy the new conditions. It will be assumed, as before, that the space 
derivatives of o, may be neglected. 

It is convenient to use cylindrical coordinates x=rcosa, y=rsing and z=2, 
the z axis being so directed that « increases by clockwise rotation when observed 
along Oz. In these axes we have assumed 


tae os 
H,= Ox =Jr=Jz=% 


so that only the « component of equation (3) need be considered. Let us define 
new components of current and electric force 7, and F,, perpendicular to H 
and lying in planes through Oz. Obviously j, is zero and F, and F, are the 
only components of electric force, so that equation (6) takes the form 


Jx=%F ,—02F,, 


0=0,F,+0,F,, 
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from which j,=03;F,. With equation (3) (neglecting displacement current) 
this gives 

OH, OH, _ , 
Oz nie “Or = 4703F,. (27) 





Multiply both sides by 7 and differentiate with regard to r: 


0 (0H, _OH, 7) 
(+ Fe =) = 41703 > (rF,) 











or 
= 4707 curl, F. 
Together with equations (10) and (14) this gives 
0/.0H, 0H,\ _ oH, 2a 
> (r <0 ) = 4703 {r — += (rv,H, - reHt,)} : (28) 


This equation (apart from the omitted displacement current term) is simply 
one of the component equations of equation (18). The second may be similarly 
derived by differentiating equation (27) with regard to z and proceeding as 


before; it is 
@ (dH, dH, oH, 2 
s(% ae) 47001 7 5 (Hl, oll). 


Thus for the much more complex magnetic field system we have 
2H = 4770, {‘F —curl(v x Hn} . (29) 


Only the meridian plane components of this equation are applicable, but these 
are sufficient for the solution of many problems. In such cases the gas behaves, 
as far as electrical induction effects are concerned, as though the gas were entirely 
unaffected by the presence of the magnetic field, its conductivity remaining op. 

A special case of the above configuration of magnetic field occurs when the 
magnetic field is everywhere parallel to the axis Oz and conditions are uniform 
along Oz. Theelectromagnetic field is sufficiently well described by equation (28), 


which simplifies to 
0 (0H 0H a 
or (- =) = 470, {r Fv + x (ro,H)}; 


where H is now written for H,,. 

These results, summarized in equation (29), will only be applied to one 
relatively simple problem: the degree to which solar and interstellar magnetic 
fields are ‘‘frozen into” the conducting gas. The degree is generally assumed 
to be very high, although the reason is not clear. In fact, if the Hall conductivity 
is neglected, the assumption is usually not justified. 

When the meridian plane component of v vanishes, the meridian plane 
component of equation (29) reduces to \/7H=470;(0H/dt). Cowling (15) has 
derived a similar equation relating to the decay of a sunspot field at photospheric 
levels where o, may be neglected, so that o,=0,=0,=0. He points out that 
the time of decay is of the order o/?, where / is the radius of the spot field. For 
a field of radius 10° cm and o= 10-8 e.m.u., the time is about 300 years, so that 
the field is effectively frozen in. If Hall conductivity were negligible in the 
chromosphere and corona (as at lower levels) the times of decay of fields of 
similar dimensions would be of the order 7 sec and 0-005 sec respectively, 
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The much shorter times are due to the greatly reduced value of o, (=o, when 
Hall current is neglected) at these higher levels. The effect of Hall current is 
to increase these times to about 10* years and 2x 10° years respectively. It 
should be noted that the increase is much greater than if o, itself replaced o, in 
equation (29). 

For interstellar HII regions (as in Table I) a field of dimensions (say) 
10° parsecs would have decay times of about 1 year and 10! years depending 
on the absence or presence of Hall conductivity. The time remains at 10! years 
even if Hall current is prevented from flowing, for the direct conductivity then 
rises from o, to 03. It is evident that the magnetic fields concerned are frozen in 
mainly because of Hall conductivity. 

The problem of axially symmetric poloidal fields with Hall conductivity 
has been considered previously by Ferraro and Unthank (16) and by Sweet (17). 
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‘THE MOTION OF IONIZED GAS IN COMBINED MAGNETIC, 
ELECTRIC AND MECHANICAL FIELDS OF FORCE 


J. H. Piddington 
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Summary 


Transient electric and mechanical forces uniform in space are applied to 
a gas in the presence of a steady magnetic field. The current transport 
equation for the anisotropically conducting medium is used to determine 
the subsequent motion of the gas, the internal electric field and the current 
density. ‘These are damped oscillatory functions leading to a steady state. 

There is a close connection between the effects of the electric and 
mechanical forces: a steady internal electric field perpendicular to the 
magnetic field cannot exist unless accompanied by a mechanical force. The 
relationship between the two is examined. Neglect of this relationship has 
already proved fatal to many astrophysical theories. On the other hand, 
there seems little justification for treating mechanical forces as “‘ equivalent 
electric fields’’, as is often done. 

This and a companion paper appear to provide, for the first time, a 
satisfactory system for using transport equations in the study of dynamical 
problems in nearly uniform, anisotropically conducting gases. Such 
problems are proving increasingly numerous and important in astrophysics. 





1. Introduction.—The motion of ionized gas in fields of force may be studied 
by summing individual ion motions or by the use of transport equations. The latter 
method has many advantages and is used here. 

The transport equation of interest is the relationship between an applied 
electric field represented by the vector E’ and the resulting current density j. 
In an isotropically conducting medium and for slow changes in E’ this is simply 
Ohm’s law. When a magnetic field H is present it takes the form: 

j= oF’ + 0,E', + 0,H * E’//H, (1) 
where E’,, and E’, are the components of E’ parallel and perpendicular to H and 
Go, , and oy are the the components of the conductivity tensor. The current 
component o,E’, is not influenced by the presence of the magnetic field and it 
provides no mechanical reaction on the gas. Without loss of generality it is 
neglected in the present paper. The other current components are the ‘‘ direct ” 
and Hall currents in that order. 

The current given by equation (1) causes a whole chain of associated events: 
first a distortion of the magnetic field according to the equation : 


curl H= 47j. (2) 

Electromagnetic units are used throughout and the displacement current neglected. 

The current also provides a mechanical reaction j x H on the gas and this, with 

hydrodynamic parameters, determines the gas velocity v. The distortion of the 

magnetic field and the gas movement in their turn help determine E’, which is the 
electric field experienced by the gas: 

E’=E+E,+u~x i, (3) 
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Here E is a potential field due to a distribution of external electric charge, E; is 
an induction field due to changes in the distribution of H and the last term is the 
electric force per unit charge due to movement of the gas. The parameters 
j, H, vand E are thus linked in a complete chain of cause and effect; for example, 
a steady electric field E’ cannot be maintained in an ionized gas unless accompanied 
by a mechanical force of appropriate magnitude and direction. 

Failure to consider every link in this chain has proved disastrous for many 
astrophysical theories, some of which have been discussed (1) already. 

Although in principle j, H, v and E’ may be determined simultaneously, 
this is often difficult or impossible. The problem may be greatly simplified, while 
reality is still retained, by neglecting the effect of j on the magnetic field. In a 
separate paper (2) this effect is discussed, but at present it is neglected so that 
equation (3) reduces to 


E’=E+vx i. (4) 


This procedure has been used by Cowling (3, 4) in studying the motion of solar 
gas across a magnetic field. The results only apply when dealing with regions 
with linear dimensions small compared with H/47j. In larger regions induction 
shielding is effective. 

The values of v, j and E’ are determined for a gas in combined magnetic, 
electric and mechanical fields of force. A transient treatment is used so that con- 
ditions prior to the steady state may be studied. The results are discussed in 
connection with different types of ionized gas which show marked and, at first 
sight, surprising differences in dynamical properties. The similarity between 
the effects of mechanical and electric forces applied to an ionized gas have led 
some workers to treat the former as “‘ equivalent electric fields ” (see reference 
(4) for example). This approach has led to some confusion and seems neither 
necessary nor desirable. 

This paper and its companion (2) appear to provide, for the first time, a com- 
plete and self-consistent system for using transport equations to study dynamical 
problems in nearly uniform, anisotropically conducting gas. Such problems are 
profuse in the field of astrophysics. Most of the groundwork was already available, 
of course, due in very large part to Cowling (3, 4, 5). The required degree of 
uniformity of the gas is discussed in the companion paper and below; departure 
from uniformity may introduce boundary problems, which will be discussed 
elsewhere. 

2. The motion of ionized gas in a steady magnetic field and transient electric and 
mechanical fields—Consider a uniform ionized gas permeated by a steady, 
homogeneous magnetic field of strength H. We confine our attention to electrical 
and mechanical fields and motions perpendicular to H. Take rectangular right- 
handed coordinates with the z axis parallel to the magnetic field. At time t=o 
apply an external (potential) electric field with components FE, and E, and alsoa 
mechanical force F dynes cm-* along Ox. ‘This constitutes the most general 
arrangement of (uniform) fields of force. 

A transient treatment is necessary to study conditions prior to the steady 
state; this study is necessary for a full understanding of the gas dynamics involved. 
It might well be asked whether there is any physical significance in the transient 
treatment. This is because of the statistical nature of o, and o, which do not 
determine the current at an instant but only over a period long enough to average 
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random-ion motions. The question is answered in Section 5 where it is shown 
that the treatment is permissible provided time t=o really represents a short 
interval of time commencing slightly after the application of the electric field. 

At time ¢ the gas velocity components are v, and wv, and the electric field 
components (as experienced by the gas) are (E,+v,H) and (E,—v,H). Hence, 
from equation (1), the current components are : 

je=o,(E,+v,H)—0,fE,—v,H), (5) 

Jy=o,(E, —v,H)+0,{E,+v,H). (6) 
Under the influence of the force F and the electromagnetic force j x H the gas, 
of mass density p, moves according to: 


pv,=jyH + F=Ho,(E,—v,H)+ Ho f{E,+v,H)+ F, 
pv, = —j,H = — Ho,(E,+v,H)+Ho{E,—v,H). 
‘These may be rearranged and solved to give: 
v,=E,/H+ F/H*o; +e (A, cos Bt + f, sin ft), (7) 
vy = —E,/H — Fo,/H*o,0, + eA, cos Bt + B, sin Bt), (8) 
where «= H*o,/p, B= Ha,/p and o,=0, + 0,2/0,, being the value assumed by the 
direct conductivity (normally o,) when Hall current is prevented from flowing (5). 
The constants A, and A, are chosen so that v,=v,=0 when t=o, so that 
A,= —E,/H— F/H*o; and A,=E,/H + Fo,/H*o,0;. The constants B, and B, 
are found by differentiating equations (7) and (8) with respect to time and com- 
paring with equations (5) and (6) at zero time. This gives B, = A, and B,= — A,, 
so that: 
v,=A(1—e™ cos Bt) + Be sin Bt, (9) 
v,= — B(1—e™ cos Bt) + Ae “sin Bt, (10) 
where A=E,/H+ F/H*o; and B=E,,/H + Fo,/H*o,03. 
The electric field (or force) experienced by the gas is given by equation (4): 
E’ 7. = — Fo,/Ho,o,;+ He*(B cos Bt + A sin Bt), (11) 
E',, = — F/Ho,+ He*(A cos Bt—B sin ft). (12) 
The currents are given by equations (5), (6), (g) and (10): 
Je=(o,E,—o2E,)e cos Bt + (o.E, +0,E, + F/H)e™' sin ft, (13) 
Jy= — F/H+(o,E,+0,E, + F/H)e “cos Bt —(o, E,—0,E,)e™' sin Bt. (14) 
The dynamic behaviour of an ionized gas in the fields of force may be illustrated 
graphically. The various parameters are represented by lightly damped oscilla- 
tions when f>« and are non-oscillatory when B<«. An intermediate con- 
dition best illustrates the effects and we choose 8/« = o,/0, = 5, so that o3/03= 5:2. 
The units used are chosen more or less arbitrarily: the time unit is such that 
B=1anda=o-2. The effect is first studied of a force F of magnitude such that 
F/H*o,=1(E,=£E, =0)suddenlyappliedalong Ox. In Fig. 1(a)thelocus of theend 
of the velocity vector v is shown moving in the xy plane in a spiral towards its steady 
state position S. In Fig. 1(b) the path traced by a small volume of the gas is plotted 
in the xy plane, the position being found by integrating equations (9) and (10). 
Both of these diagrams apply equally well for an electric field which gives the same 
values of A and B. That is, for EF, = Fo,/Ho,o;=5H and E, = F/H*e,=H. This 
equivalent electric field is called E, and is shown as a vector in Figs. 1(a) and 1(). 
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to the force F (along Ox) ; these are plotted against time and denoted E£’ ,( F) 


and E’,(F). The corresponding curves for the equivalent electric field E, are 


In Fig. 2(@) are shown the internal electric field components E’, and E’, due 
denoted by E’ (Ey) and E’{E,). Their steady state values are zero. In Fig. 2(d) 


the two sets of current density components are plotted similarly, being denoted 
by j.(F), j,(F) and j,(E,), j,(Ey) for the applied fields F and E, respectively. 
The steady states of all but j,(F’) are zero. 
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(b) 
Fic. 1.—The movement of ionized gas under the influence of a steady magnetic field 
(vertically outward from the paper) and either a suddenly applied mechanical force F or a suddenly 
applied electric field Ep as shown. 


(a) The locus of the end of the velocity vector v. (b) The path traced by a small volume of the gas. 


v and j throughout the medium. 


It has been assumed above that conditions are everywhere uniform. If this is 
not the case, then attention must be paid to the possible effects of variation of 


In some problems such variation is not im- 
portant; in others, however, it results in accumulations of one or more types 
of particle (charged or uncharged) in certain regions. This is a boundary problem 
and may add greatly to the difficulty of solution. 


considered in relation to each individual problem. 


For example, if the con- 


Boundary effects must be 
ductivity changes at a plane, there will be a tendency for space charge to collect 
on that plane and the values of E,, and E,, in equations (5) to (14) would be functions 
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(b) 


Fic. 2.—(a) Internal electric field components E’z(F) and E’,(F) due to a suddenly applied 
force F and internal electric field components E’ (Ep) and E’y(Ep) due to the equivalent electric 
field Ep. All plotted against time. (b) Electric current density components jz(F) and j,(F) 
due to a suddenly applied force F and electric current density components jx(Ep) and jy(Ep) 
due to the equivalent electric field Ep. All plotted against time. 
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of j, andj, and of time. Mention is made below of pressure gradient as one type 
of force acting on the element of gas concerned. It is tacitly assumed that 
boundary effects are not introduced by the gradient. 

Further discussion of equations (9)-(14) in connection with specific types of 
ionized gas is given in Section 5. 

A special case of the above problem occurs when the external electric and 
mechanical fields are zero but the gas has an initial velocity, say v along om. 
Equations (5) and (6) still hold but E,=E,=F=o and also when t=o 
U,=U, ly=0, ¥,= —0,H*v/p and v,= — oglT*vp. From these boundary con- 
ditions A,= = ~ee and hahas, so that: 


v,=ve'cos Bt and v, = — ve sin Bt. 


This result was obtained by Cowling (5) who was studying the motion of solar 
material across magnetic fields. He concluded that the motion decays in a time 
comparable with «~ or p/o, H®. Cowling was concerned with gas near the reversing 
layer at which level o, is of the same order as or smaller than o,, so that B <«. 
The motion of the gas is then too highly damped to be oscillatory. It may be 
noted that at higher levels where o,>0, the motion is oscillatory, the forward 
motion being substantially changed in time comparable with B- or p/o.H? 
which is much smaller than the time p/o,H? taken for all motion to disappear. 

3. Steady-state conditions —The transient treatment used above indicates 
oscillatory movement of the gas with angular frequency B= H?a./p. It is evident 
that if the fields of force are applied gradually over a period of about 27/8 or more, 
there will be no oscillatory motion and the steady state is assumed immediately. 
In any case this state is reached after time of the order «1, owing to the ex- 
ponential decay term. 

The steady-state conditions are: 


v,=E,/H+ F/H’o3, v,= —E,/H— Fo,/H*o,03, (15) 
E,'= — Fo,/Ho,0;,  E',= — F/ Hos, (16) 


Je=0, jy=—F/H. (17) 

When the force F is zero the gas moves with velocity E/H perpendicular to 
the applied electric field E. It experiences no internal electric field, no current 
and no energy dissipation. 

When the externalelectric field is zero the gas moves with velocity F/{H?(o,03)"?} 
in a direction ¢, measured clockwise from F, where tan ¢=o,/0,. At the same 
time an internal (induction) electric field is experienced by the gas of magnitude 
E’ = F/{H(o,0,)"*} at an angle 7/2+¢ measured clockwise from F. ‘The com- 
ponents and vectors of velocity and internal electric field are shown in Fig. 3. 
which is drawn for a gas having o, somewhat larger than o,. 

In the case of a very tenuous gas when collisions may almost be neglected 
o,/0,;>1, so that d&7/2. The gas moves almost perpendicularly to the force 
which causes its motion and E’ is nearly antiparallel to F. 

It is important to observe the close connection between the internal electric 
field and the mechanical force: one cannot exist without the other. If an external 
electric field is applied, the gas moves so as to neutralize this by an induction field. 
In so doing it may experience a force due, say, to the development of a pressure 
gradient. The steady condition is shown in Fig. 3 by the vectors F and E£’. If 
no steady force opposes gas motion, then the internal electric field dwindles to 
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zero. Again, if a mechanical force is applied without an electric field the gas 
moves until an induction field provides a balance as in Fig. 3. Failure to observe 
these interdependent factors has led to much confusion in the past. 

To complete Fig. 3 an electric field vector Ey, is drawn equal and opposite to 
E’ and so at an angle tan-! o,/c, to F. This electric field is exactly equivalent to 
F in determining gas velocity, both transient and steady state. They are not, 
however, equivalent in other respects and should not be treated as such. Thus 
from equations (16) and (17) and Figs. 1 and 2 they are clearly not equivalent in 
determining the values of internal electric field and current. 

Steady-state conditions for the case of a lightly ionized gas have been discussed 
by Martyn (6) who has determined the ion and electron drift-velocity vectors 
due to an applied electric field. He finds an internal force due to collisions 
between the ion-electric gas and the neutral gas. ‘This force bears the same 
relationship to the electric vector as that given in Fig. 3 (when allowance is made 
for the different assumed directions of the magnetic field). Following the pro- 
cedure used here of considering the gas as a whole (neutral as well as charged 
particles), an external force of the same magnitude must be exerted on the neutral 
particles to prevent them from moving. This corresponds to F in Fig. 3. 





Fic. 3.—Vector diagram of steady-state gas velocity v, internal electric field E’ and current 
j for a gas in a magnetic field (vertically outward from the paper) and subject to a mechanical 
force Fcm~-* as shown. The electric field Ep shown causes the same movement of gas but no 
internal electric field or current. 


The rate at which the force F does work on the gas is Fv, or F(E,/H + F/H*o;). 
The second term corresponds to the Joule heat dissipation F?/H*o; =j?/0,;=jE’, 
and shows that o, is the conductivity which determines rate of energy dissipation, 
an important conclusion reached elsewhere (2). The first term, FE,/H = —j,E,, 
corresponds to the energy supplied by the force F and transmitted by the gas to 
the external electrical system responsible for supplying the electric field E,. 
‘This energy may result in a build-up of electric polarization at boundaries or in 
the dissipation of Joule heat in external circuits. 

A case which may have some astrophysical interest is that of a gas moving 
with velocity v,=£,/H under the influence of E,. If this velocity is destroyed 
by a sufficiently powerful force antiparallel to Ox, the gas develops a velocity 
vy =(e2/0,)v,. The value of 04/0, is often large, being 10°-10° in the Sun’s 
atmosphere and in interstellar H 11 regions, so that even relatively small values 
of v, and therefore of E, may result in values of v, approaching the velocity of 
light. In practical problems close consideration must be given to the build-up 
of electric polarization at boundaries. This effect tends to prevent the gas mov- 
ing as above. 
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Before proceeding with, and in order to facilitate the discussion of the equations 
of Section 2, expressions for the conductivity coefficients of some gases will be 
introduced. 

4. Ionized gas conductivity.—The conductivity coefficients of equation (1) 
have been determined for various gas mixtures by Chapman and Cowling (7) and 
Cowling (3). Itis proposed to discuss two extreme types: a binary gas consisting 
of electrons and one type of heavy ion, and a lightly ionized gas containing the 
same two ingredients with a much larger number of neutral atoms or molecules. 

In a binary gas: 


Gg=n,e?7/m, 0, =09/(1 +w*7?), o,=agw7/(I+wr?) and o;=a9 


in electromagnetic units, where n,, e, m and 7 are the electron density, charge, 
mass and mean interval between collisions with ions and w=He/m. It was first 
shown by Cowling (§) that o,==0, +0,”/o, was the value assumed by the direct 
conductivity (normally o,) when the Hall current was prevented from flowing 
by electric polarization. Although Cowling was discussing a binary gas there is 
nothing restrictive in his argument which applies for any gas. ‘The importance 
of a; in some ionospheric problems has been shown by Baker and Martyn (8) 
and others. 

The important identity o)=o, holds only for binary and other gases in which 
collisions by a single component (electrons) with one or more other components 
determines the conductivity. For such gases r oc n, 17%", where T is the gas 
temperature, so that o, and o, depend only on 7' and not on”, or H. ‘This greatly 
simplifies many problems as already shown (2). 

The conductivity components of a lightly ionized gas each have two terms 
corresponding to collisions between electrons and neutral atoms and between 
ions and neutral atoms. Only one type of ion is assumed and collisions between 
electrons and ions are neglected. We have: 


Oy =0,+0;=n,e?7,/m+n; Ze 7;/M, 
oy =9,|(t +072) +0;/(1 + 2272), 
02 = 0,7, /(1 + w*7,”) —0,0,7;/(1 + 2"7,7), 


where 7, and 7; are the collision periods of electrons and ions with the neutral 
molecules and n;, Ze and Q are the ion density, charge and gyro-frequency 
respectively. 

A difference of fundamental importance occurs in the values of og and o; for 
the above gases when they are very tenuous or the magnetic field very strong, 
that is when w7, wr, and 7; are all large. For the binary gas o,, o, and o; may 
be written ,e/H .1/w7,n,e/Handn,e/H.wr, so that o3/0,=0,/0,=w7 which is 
large compared to unity. For the lightly ionized gas, which is neutral so that 
n,e=n,Ze, we have: 


ar + iam and o?=o0, so that 
o1=— da“=0, s at o; =o}. 
I H\ wr, | 7; : — 


‘Thus G3/0,>0 and oya,>0 and a3/49=0,/09= 1/(w7£27;), 


which is very small instead of unity for the binary gas. ‘This profound, and 
at first sight surprising, difference in behaviour of the two gases is discussed 
in Section 6. 
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Because of their great importance in astrophysics “ impure binary gases ” 
also require discussion. One such gas comprises electrons and ions together 
with some neutral particles (called atoms for brevity), but in which the electrons 
collide much more frequently with ions than with atoms. Because of the smaller 
collision cross-section of atoms they may be much more numerous than ions 
while still satisfying this condition. ‘The conductivity components of such a gas 
are the same as for binary gas provided one vital condition is satisfied; this is 
discussed in the following section. 

Another gas which might be called impure binary comprises electrons and 
two or more types of ion, in which the contribution to conductivity due to collisions 
between ions is negligible. This gas has the same conductivity as a binary gas. 

5. Discussion of equations (9)(14).—The quantities o, and o, determine the 
average current which flows when an electric field is applied to a medium for a 
sufficiently long period. Generally this period is taken as long compared with 
both the gyro-period and the collision period. In Section 2 this rule was not 
observed and the physical significance of the treatment needs discussion. 

5.1. Binary gas.—Consider first a binary gas in which wr>1; the con- 
ductivity is provided by the electrons and the correct values of o, and a, are 
established in a period during which the electrons have made several gyrations. 
In this period the ions will hardly have moved and may be considered at rest. 
Subsequent movement of the ions constitutes movement of the gas as a whole 
and introduces the term v x H into the equations, thus changing j but not o, 
or gg. It is not necessary that the period during which j is averaged should be 
long compared with the collision period. Each collision by any electron con- 
tributes to j and it is only necessary that the total number of collisions of all the 
electrons in a volume small compared with the total volume of the gas should be 
statistically large. Hence when wr>1, o, and o, determine the average current 
flowing in a period short compared with the gyro-period of the ions. 

A better understanding of the result above may be obtained by assuming a 
force F applied to the gas simultaneously with the electric field. Let F be a 
gravitational force acting almost entirely on the ions and of such strength that 
it exactly neutralizes the electric force due to the applied electric field. Neglecting 
thermal motion the ions remain stationary while the electrons move under the 
influence of the electric field. Stable conditions are established within a few 
electron gyro-periods, the electrons conducting current of density j such that 
jx H=F. When F is not present the conditions are only changed in so far as 
F is replaced by the force necessary to overcome inertia of the protons. 

When the collision period is not long compared with the electron gyro- 
period, the necessary averaging period is not increased although it may be 
decreased if wr<1. Itis concluded that for a binary gas the treatment of Section 2 
is valid provided “‘ zero time” really means a time interval during which the 
electrons may perform several gyrations. 

The binary gas equations show that o2/o,=w7, so that 8 =aw7, and when wr 
is substantially greater than unity equations (11)-(16) represent an oscillatory 
condition. ‘The angular frequency § is found from the conductivity equations 
to be equal to the angular frequency of rotation of the heavy ions in the magnetic 
field, HZe/M. ‘The amplitude of the oscillation is found for the case where 
F=E,=0, so that 4=E,/Hand B=o. Since we have shown the correspondence 
between a mechanical force and an electric field, this example is sufficiently general. 
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The amplitude of oscillation (found by integrating equations (g) and (10)) is 
A/B=E,M/H*Ze along both x and y axes. There is also a steady drift with 
velocity E,/H along Ox. ‘The oscillations are damped by a factor exp (—1) in 
time «| =p/H?o,. 

When the collision period is equal to or smaller than the electron gyro-period 
B <a and no oscillations occur. Steady-state motion is acquired, as before, in 
time of the order «1. 

5.2. Impure binary gas.—A binary gas to which some atoms have been added 
may be considered as a pure binary gas moving freely within an atomic gas. All 
the considerations above regarding a binary gas are equally applicable, the value 
of p corresponding to the binary gas only and not the whole gas, and F being the 
force acting on the binary gas. However, when the steady state given by equations 
(15)-(17) has been attained, a further process takes place. The collisions between 
the binary and the atomic gas, although few, must sometimes be taken into 
account. This is done by calculating the rate of transfer of momentum from the 
binary to the atomic gas, namely (”,m/7, +n,M/7;) (v— a), where v and uw are the 
velocities of the binary and atomic gases and 7, and 7; are defined as for a lightly 
ionized gas. Although the electrons and ions have drift components other than 
v, these components make no contribution to momentum transfer. This has 
been demonstrated by Martyn (6) for a lightly ionized gas and is equally applicable 
here. ‘The importance of the frictional force between the gases depends on the 
relative values of 7, and 7; and the time during which the fields of force are applied. 
Thus for oscillatory fields of period <7; it may be neglected. If the period 
is comparable with or greater than 7, , it may be important. 

Gas containing two types of heavy ion but having no appreciable ion 
conductivity is really a mixture of two binary gases. In the oscillatory stage each 
gas moves approximately independently with the gyro-frequency of its heavy 
ions. In the steady state the gas may be considered as a single binary gas with 
appropriate mean mass density. 

5.3. Lightly ionized gas.—In a lightly ionized gas, electrons sometimes 
contribute nearly all of the conductivity, in which case the treatment used in 
Section 2 is valid. If, on the other hand, ions contribute appreciably, then the 
treatment is not valid. In this case ion transient movement may be studied by a 
consideration of the electrons and ions as constituting a binary gas with infinite 
collision period, moving within an atomic gas. The conductivity components 
of the binary gas are: 


o,=0, o,=1,e/H, o3=09= ©. 


‘lhe constants 4 and B of equations (g) to (13) become E,,/ H and (E.,/H + F' Ha.) 
anda=o. If wr, and Q7, is both much greater than unity then the binary gas 
oscillates according to the equations, F being the force on the binary gas only. 

If wr, or Q7,; is comparable with or smaller than unity the binary gas 
does not oscillate. In any case after a time of order wr,/Q or 7;, whichever is the 
smaller, the steady state is attained. ‘The motion is then described by equations 
(15), the electric field experienced by the binary gas and the current flowing are 
given by equations (16) and (17). The force F is the frictional force 
—(v—u)(n,m/r,+n;,M 7;) between the binary and atomic gases. The latter moves 
under the influence of a force — F together with any other forces acting directly on 
the atoms. 
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For steady-state conditions a lightly ionized gas is best considered as a whole, 
relative internal movement of the binary gas being ignored. The value of v then 
refers very nearly to the atomic gas (this is discussed further in the next section) 
and E’ and F also refer to this gas. ‘The conductivity components are those of a 
lightly ionized and not a binary gas. 

6. Dynamics of ionized gas.—There is as yet no complete and self-consistent 
system of using the conductivity equation to solve gas dynamical problems. 
The reason is that different workers have developed different methods of approach 
which are difficult to reconcile with one another and in some cases are faulty. 
Methods tend to differ when dealing with different kinds of gas or different 
combinations of mechanical and electrical forces. Three aspects of the problem 
require discussion in which the present approach is compared with some 
others. 

6.1. The gas velocity.—A point which has caused some confusion is the 
definition of v in the equation j=[c](E+ux H), [oc] being the conductivity 
tensor and E the electric force experienced by unit charge at rest. The concept 
of conductivity was introduced as a property of solids and it was evident that v 
should refer to the solid itself, even though charged particles moving within the 
solid had different velocities and so, individually, experienced different electric 
forces. In the case of gases and, in particular, of a fully ionized gas, the position 
is not so clear because then all the particles are electrically charged and each type 
of particle moves with a different average velocity. However, the value of v, 
as pointed out by Cowling (4) is the “‘ velocity of the mass as a whole”. It might 
help to be even more specific and define v as the net momentum per unit mass 
of the medium. The different drift velocities of the components are cared for 
in the derivation of the conductivity. 

It is the definition of v which accounts for the extreme difference in [o] of 
very tenuous binary and lightly ionized gases mentioned in Section 4. Students 
of lightly ionized gases sometimes tend to regard binary gases as extreme cases of 
lightly ionized gases, from which most of the neutral atoms have been removed. 
This is not so, because of the different definition of v for the two gases. In the 
lightly ionized gases v is (very nearly) the velocity of the neutral atoms, while in 
a binary gas it is (very nearly) the velocity of the heavy ions. Thus in a lightly 
ionized gas, in which the atoms remain at rest, movement of ions modifies j 
without changing v—hence it must affect the value of [oc]. In a binary gas, on the 
other hand, a velocity v of the ions corresponds almost exactly to a velocity v of 
the whole medium and thus determines the electric force experienced by the 
medium but not its conductivity. The apparent anomaly is now easily explained. 
When an electric field E, is applied perpendicular to H in a lightly ionized gas in 
which collisions may be neglected, the ions and electrons move separately as a 
binary gas and assume a velocity such that E, + v x His zero andj iszero. Since 
the electric field E, experienced by the gas as a whole (the neutral atoms) is not 
diminished it follows that o, and oc, must be zero. In the case of a binary gas, when 
E, is applied without a mechanical force, the electric force E, + v x H experienced 
by the whole gas falls to zero. Hence j is zero, not because o, and a, are zero 
but because the total electric force is zero. If a mechanical force restrains the 
ions from moving, a (Hall) current flows, indicating a finite value of op. 

6.2. Pressure gradients.—The results of Sections 2 and 3 are derived for a 
uniform gas. However, provided conditions do not change appreciably in the 
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region considered, they may be applied where pressure and density gradients 
exist. ‘he various parameters are then given by the equations of Sections 2 and 
3 in which F may now contain a term in grad p, where p is the total gas pressure. 

If F includes other forces, say gravitational, these may be lumped together with 
grad p. It may seem surprising that a gravitational force which acts almost ex- 
clusively on the heavy ions (in the case of a fully ionized gas) is equivalent to a 
pressure gradient which acts on both ions and electrons. The reason is that a 
space-charge electric field is set up which effectively links the two gases. If there 
is a pressure gradient there is, in general, a density gradient, and as the electrons 
have less inertia they are pushed ahead of the heavy ions in the direction of 
decreasing density. Only a small movement is generally sufficient to create a 
space-charge field sufficient to transfer the necessary momentum from the 
electron to the ion gas. The space-charge electric field may be considered as an 
internal part of the whole mechanism and need not appear in the gas-dynamical 
equations unless it is of interest in its own right. 

In the case of a partially ionized gas there may be very different forces acting 
not only on electrons and heavy ions but also on neutral particles. Thus radiation 
pressure might act exclusively on the latter while electric forces acted only on the 
charged particles. ‘This would seem to necessitate great complexity in the 
gas-dynamical equations. However, only in the transient stage need the separate 
motions of the gas components be considered. In the steady state an appropriate 
drift of neutral particles relative to ions provides a frictional force between them 
which allows a suitable transfer of momentum as in the case of the space-charge 
field. The relative motion and frictional force may remain an internal matter 
not considered in dynamical problems. 

It is seen that whenever possible any ionized gas should be considered 
dynamically as a single gas, relative motions of its components being ignored as 
well as internal electrostatic fields. In the case of a partially ionized gas this can 
only be done provided changes are slow enough'to allow a steady drift to be 
attained. 

6.3. Equivalent electric fields—Cowling (4), Schliiter (9) and others have 
developed gas-dynamical equations from the individual equations of motion of 
the electron gas or the positive ion gas or both. Cowling finds the total force on 
the electron gas to be 

grad p,—n,e[E+(v+V) x H]+n,mV, 

where p, and V are the partial pressure of the electron gas and its velocity relative 
to the heavy ion gas (a binary gas is being considered), v is the collision frequency 
of the electrons, and the other terms have been defined above. He then deduces 
a total ‘‘ effective ” electric intensity on the moving gas 

E’=E+vux H+ grad p,/n,e. (18) 
If E’ is perpendicular to H, then it is the field which determines the total current 
density according to the equation: 


, 


j=0,E' +oHx%, (19) 


Thus the “ effective ’ electric intensity includes a force term grad p,/n which is 
termed an “‘ equivalent ”’ electric field. 

In subsection 6.2 above it was stated that a pressure gradient could be lumped 
with a gravitational (or any other) force in determining the gas motion and electric 
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current density. In Section 3 it was shown that a force F caused an internal 
electric field of magnitude F/H(c,0;)"* which would be effective in causing 
current flow according to equation (19). For low values of collision frequency 
in a binary gas: 
H(c,03)'?=n,e, 

so that the “ equivalent” electric field is approximately F/n,e or grad p/n,e 
when F comprises a pressure gradient. ‘This appears inconsistent with Cowling’s 
result in which only the partial pressure of the electron gas is introduced. 

To resolve this inconsistency it is necessary to consider the actual physical 
effects of the pressure gradient. ‘That part of the pressure gradient which acts on 
the electron gas causes this gas to move slightly away from the heavy ion gas 
and so creates a space-charge field which may be included in the term E of 
equation (18) or kept separate and called an equivalent electric field. ‘This 
electric field is present when the steady magnetic field is zero. ‘That part of the 
pressure gradient which acts on the heavy ions may be considered as acting on the 
gas as a whole (as would a gravitational force), dragging it across the magnetic 
field. ‘The motion induces an electric field which may be included in the term 
v x H or kept separate and called an equivalent electric field. This component 
occurs only when there is a steady magnetic field with a component perpen- 
dicular to grad p. ‘The two electric fields due to the pressure gradient may 
be kept separate if desired, or combined to give a total equivalent electric field 
grad p/n, e. 

Schliiter (9) also found an “‘ equivalent’ electric field which differed from 
that of Cowling and from the value given above. 

The concept of ‘‘ equivalent ”’ electric forces may be useful but it may also 
introduce some confusion. ‘The alternative approach used above is to deal only 
with real electric fields, all of which are included in the terms of equation (4), 
and with mechanical forces which help determine these fields as well as the values 
of gas velocity and current density. 
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Summary 


The possibility is examined of a pair of originally single stars becoming 
gravitationally bound together as a result of their interaction with interstellar 
material. In the first instance, a special model is considered in which it is 
assumed that this interaction is in accordance with the Bondi—Hoyle- 
Lyttleton accretion theory. It is found that in this model such binary star 
formations are possible, but under certain conditions which could have 
existed only at an early stage in the evolution of the stars. In the light of the 
results of this preliminary examination, a more realistic model is considered 
and this leads to a very effective process for binary formation when the same 
properties of the interstellar material are assumed as in the accretion theory. 





1. Introduction.—There appear to be only three logically possible theories 
of the origin of binary stars: the binary systems may have been formed by the 
disintegration of originally single stars, they may have been born or created as 
double stars, or they may have been formed by the coming together of two 
originally single stars. 

We shall consider briefly the three types of theory. The first or fission theory 
was for some time the accepted one. However, later examination (1, 2) indicated 
that it is inadequate to explain the formation of binaries. According to the second 
theory, which is the one now receiving most attention, binary stars were formed 
by the condensation of two stars within each other’s gravitational influence. 
In the third theory, two originally independent stars are considered to come 
together and remain together. It can be shown that if two stars approach from a 
great distance, they must necessarily separate again to a great distance in the 
absence of forces other than their gravitational attraction. The most obvious way 
to introduce another force is by the presence of a third star. ‘Thus, if three stars 
approach from a great distance, it is possible for one to recede to infinity leaving 
the other two gravitationally bound together. 

The first theory we wish to consider here is of the third type but the place 
of the third star is taken by interstellar material. In other words, two stars are 
considered to approach from a great distance apart in a region containing inter- 
stellar material. During the approach and subsequent encounter, each star will 
experience a resistance to its motion and so will lose energy. If this energy loss 
is sufficient, the stars will not again recede to a great distance apart but will remain 
gravitationally bound together. Once the binary has formed in this way, its 
orbital elements may be subsequently affected by the interstellar material, but it is 
not proposed to consider this in the present paper, 
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If we examine more closely the second type of theory of binary formation, 
we see that some similar process must occur, because after the stars have con- 
densed their subsequent motions must for some time be influenced by the star- 
forming medium. It was during this early stage that we consider the mechanisms 
here put forward to have been of most importance. 

In Section 2 of this paper an elementary mathematical treatment of such a 
binary-forming encounter is given. This treatment is based on the assumption 
that the stars experience a mass increase in accordance with the Bondi-Hoyle- 
Lyttleton accretion theory. In Section 3 the result of this treatment is discussed 
and its implications are considered. In Section 4 a limitation of the elementary 
treatment is considered and a different model is suggested that can be treated 
more realistically. ‘This leads to a mechanism which appears to be an efficient 
agency for the formation of binary stars. 

2. Elementary mathematical treatment.—In considering the dynamics of a 
binary-forming encounter between two stars, as outlined above, we must imagine 
the stars to approach from a separation of the order of the local mean interstellar 
distance, since at greater separations each star is more likely to be influenced by 
other neighbouring stars. The maximum initial velocities of the stars necessary 
to produce a capture will depend on the initial directions of motion and on the 
masses of the stars. A full investigation of the binary-forming encounter would be 
quite difficult in view of the numerous parameters involved. However, a simple 
example has been investigated in order to find the order of magnitude of the initial 
velocities. ‘The case considered is that of two stars of equal mass approaching 
with equal speeds in almost opposite directions. The word “ almost” in the 
last sentence is intended to mean that from the mathematical point of view the 
motion of the star can be considered rectilinear, but on the other hand the stars 
must miss a collision at the point of minimum separation. The centre of mass 
of the stars is taken to be at rest in the medium which has density p. Starting 
with an initial separation d, the initial speeds are found, so that stars never again 
separate to a distance of more than $d. The choice of this maximum subsequent 
separation is arbitrary but it must be such that the probability is low that the 
stars will be torn apart again by gravitational interaction with neighbouring stars. 

We must next consider the interaction between the interstellar material and 
the stars. It is proposed to assume the mechanism which has been developed in 
the accretion theory of Bondi, Hoyle and Lyttleton. For details of this, 
see for example (3) and the references therein. For our purpose, however, 
we only need to know that the theory implies that each star ‘“ accretes” or 
captures the interstellar material lying within a distance o from the star as it 
moves through the material. As shown in (4), 


o = 2pu/v?, (1) 


where » = mass of star x constant of gravitation and v=velocity of star relative 
to the interstellar material. For very small velocities it is seen that the value of o 
becomes very large. There is however a limit set by the presence of neighbouring 
stars. These limit the capture radius o to about one-half of the distance of the 
star from its nearest neighbour, because it is only within this distance that a star 
can be considered to have a predominating gravitational effect on the interstellar 
material. We are going to suppose that the velocities are so small that this maxi- 
mum capture radius applies. The rate of capture of material by the star is 


45 
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therefore 7o%pv, where a is half the distance of the star from its nearest neighbour, 
which in our case will be the other star partaking in the encounter. We must 
remember however that this capture rate applies subject to the condition that 
the velocity does not exceed the value given by (1), i.e. v<(2y:/0)"?. 

Let a frame of reference be placed such that the origin O is at the mass centre 
of the stars and let the axis Ox pass through the starting points of the stars. Then 
the motion of either star will be almost rectilinear and coincident with Ox. Let x 
be the position, v the velocity along Ox and m the mass at time ¢ of the star which 
started on the negative part of the axis Ox. Then 


dx 
a =* (2) 
The equation of motion is 
d(mv) Gm? 
tine ge (3) 


when the stars are approaching, where G is the constant of gravitation. The 
rate of accretion is 


om = nxtp|v |. (4) 


The velocity of the star under consideration is positive, so dividing (4) by (2) gives 


—— =7px’, 


dx 
whence m= 4rpx* + M (5) 


where M is the mass of the star when it reaches the origin. From (3), 


dv, dm _ Gm* 
wa Tek” Boe 
and using (2), (4) and (5), this becomes 
du anpx*v® = G(M + }7px*) 
dx + + }npx® 2x2 ; 
of which the solution is 
I M? zpM*x? = 2*p?Mx>_—s wp 3x8 
= cartel O(- ae tt ae tgs 4] © 
A being a constant of integration. Now (5) implies that as the stars approach 
each other, they accrete the material lying within the cone whose vertex is at O 
and whose semi-angle is 7/4. The stars narrowly miss a collision at O and then 
proceed and come to rest at x= +}d. Since after passing O each star will be 
entering a region which has just been cleared of material by the other star, we 
assume that the stars do not accrete between O and their first position of rest. 
Thus, in this part of the motion, we have from elementary considerations (or 
else by putting p=o and changing the sign of G in (6)), 











= —— +B, (7) 


B being a constant which is determined by the fact that v=o atx=}d. So 
=—2GM/d. The condition at the origin is that at x=5 (where 3 is a small 
positive quantity), v? given by (7) must be the same as v* at x= —5 given by (6). 
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From this it follows that A= —2GM2?/d. To obtain the maximum initial velocity 
of a star we put x= —}din(6). We obtain 
Gpd? 7 Ci 7 
= Seana carci teat A eta 
er ~zyl ait 960" * aa: | (8) 
a = 


2 73 


where « = M/pd°, which determines v for various values of M. It will be noted 
that when « is large enough, v? will be negative owing to the negative coefficient 
of «* in (8). Thus there is a limit above which « must not go. This limit is 
obtained by equating to zero the cubic expression in (8). The solution is approxi- 
mately «=0-135. The physical significance is that stars with « >0-135 cannot 
be brought to rest at a separation of }d even if their initial velocities (at separation 
d) are zero. 
TaBLe I 


Masses of stars resulting from the binary-forming encounter of stars of initial separation 
4 parsec and of given initial masses in a medium of given density 


p Initial masses (in solar masses) 
g/cm*® 0:00069 0°0069 0069 0°69 





107% 0°023* (0029) (o°09g1) (0°712) 
1077! 0°222 0°230* (0°290) (0°91) 
10-20 2°21 2°22 2°30* (2:90) 
so" 22°1 22°1 22°2 23°0* 





TABLE II 


Maximum initial velocities (km/sec) of stars in the binary-forming encounter (initial 
separation 4 parsec) 


p Initial masses (in solar masses) 
g/cm* 0°00069 600069 0069 0°69 





107% 0°0299* 
(0-0035) 


t0~™ 2°13 0°0945* 
(00035) (o-or1) 


10720 70°! 6°75 0°299* 
(00035) (oor) (07035) 


10~}9 2240 222 21°3 0°945* 
(00035) (o-or1) (07035) (0-11) 





3. Discussion of the results.—For stated values of the initial masses and of the 
density of the interstellar material, Tables I and II give the maximum initial 
velocities and resulting masses M of the stars undergoing such a binary-forming 
encounter. The initial masses are such that the diagonal figures in the tables 
(marked with an asterisk) correspond to «=0-135. For these tables the value 
of d was taken to be half a parsec. This value was taken because although the 
present mean interstellar distance in the neighbourhood of the Sun is about one 
parsec, if the binary stars were formed by the method suggested then the mean 
interstellar distance of the original single stars must have been rather less. 


45* 
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The initial velocities given in Table II are unfortunately too large for the 
assumption about the accretion rate to be satisfied. The maximum velocity for 
this assumption to hold is given by the condition (1), Applying this condition 
at the point x= —}d (i.e. putting o=4d in (1)) we obtain the velocities given in 
brackets in the table. We must therefore look upon these as representing the 
order of magnitude of the maximum velocities. It may be thought that by using 
the appropriate accretion formula for higher velocities, a higher maximum could 
be obtained. But it is unlikely that the maximum velocities could be appreciably 
increased in this way because as the velocity increases beyond the value determined 
by (1), the accretion rate falls off very rapidly. 

The bracketed velocities are very small and consequently this binary-forming 
process can only be considered to apply to stars which are almost at rest initially. 
This is one reason for thinking that this process could only have been of importance 
at the stage when the stars were condensing from the original medium. Another 
reason is the fact shown at the end of the last section, that the initial masses of 
the stars in the encounter have an upper limit. Another point arising from the 
fact that the initial velocities must be small is that in all cases the binary-forming 
encounter must be an almost rectilinear motion, since the small initial velocities 
will soon become negligible compared with the velocities developed by the 
attraction of the two stars. It follows that the treatment given in Section 2 is 
not such a special case as at first it appeared to be. 

The next point to note is that the resulting masses of the stars are almost in- 
dependent of the initial masses. This is because the mass of material accreted by 
each star is much greater than its original mass. The resulting masses are therefore 
only dependent on the density p and also on d, being proportional to d*. To obtain 
a binary consisting of two components of solar mass we would require a density 
of between 10-*! and 10°” g/cm* (with d=} parsec). This is high compared with 
the densities of most clouds observed at present, for which p does not exceed 
10°"? g/cm*, but it must be remembered that when the binaries were formed, 
considerably different conditions of the interstellar material must have prevailed. 
Of course, for a given value of p, any desired resulting mass can be obtained by a 
suitable choice of d. The importance of the present calculations is that the masses 
of observable binaries can be obtained with quite reasonable values of p and d. 

In Table I the bracketed figures are included for the cases which do not result 
in captures. It is seen, of course, that in all cases the masses of the stars are 
increased. It follows that if a star does not form a binary in its first encounter 
with another star, then it can never again do so by this process in a medium of 
the same density because if we re-enter the table at an initial mass equal to the 
new mass of the star, we shall be even further removed from the region of the 
table where binary-formation is possible. It also follows that if a binary forms, it 
cannot subsequently pick up a third star to form a tertiary system under the 
influence of the medium alone. It is however possible for three or more initial 
condensations to become mutually attracted and form a multiple system. To 
investigate the case of three condensations, we may consider three equal con- 
densations at the corners of an equilateral triangle and examine their motions 
as they approach their common centre of mass through the medium. If we do 
this, we obtain the expression (8), except that the G is increased by a factor 4/ 1/3. 
Thus, the initial velocities of approach v, as given by (8), are no smaller than for a 
binary-forming encounter. Similarly for four condensations originally at the 
corners of a regular tetrahedron. 
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3.1. Amount of interstellar material remaining —When a binary-forming 
encounter of the type discussed above occurs, the mass of material accreted by 
the two stars is about one-half of that originally contained within a sphere of 
diameter d. Consequently there will still be sufficient interstellar material left 
to form a resistive medium for further accretion effects to occur. The formation 
of multiple systems leaves somewhat less unaccreted material. 

3.2. Time required for the encounter.—In any particular case (i.e. for given initial 
masses and velocities of the stars), the time required for the stars to reach the 
origin is 


er 


-aj2 0’ 


where v is given by (6). This integral was evaluated numerically for the case of 
initial masses corresponding to «=0-135 and for zero initial velocities. This 
integration was performed with the help of the initial approximation for the 
velocity 


v? =2Gm,(x + $d)/d*, 


where mz, is the initial mass of the star, which from (5) is pd*(«—7/24). The 
resulting time was found to be g:12(Gp) 1”. For p= 10~* g/cm? this is 1-12 x 108 
years and for p=10~* g/cm? it is 1-12 x 10’ years. 

3.3. Effect of unequal initial masses.—If the masses of the condensations are 
not equal initially, let the masses at any instant be m, mg, the positions x,, x, and 
the inward directed velocities relative to the interstellar material v,, v. respectively. 
Then (3) becomes 

dmv) _ Gmym, _ d(myv) 
a § (%-z,y a@&’ 





from which it follows that at any instant 
m,V, = Mv, + Constant. (9) 


If, in this case, we consider the stars to tunnel out regions of radii equal to their 
distances from the neutral point of the gravitational field (i.e. the point where the 
gravitational attractions of m, and m, are equal and opposite) then the accretion 


rate for m, is given by 
1/2 
dm, = pv m, |x; —x,| 2 (10) 
dt *L m}2+m,'? | * 


If the condensations start from rest, the constant in (g) is zero and so (10) and the 
corresponding expression for dm,/dt are equal. Thus, during the encounter, both 
masses accrete equally and since the total mass of the accreted material greatly 
exceeds the masses of the initial condensations, it follows that the masses of the 
stars after the encounter are practically equal. If, however, the initial velocities 
are not zero, the accretion rates will be unequal and so some difference may 
be expected in the resulting masses. It is not possible to estimate the extent of 
this difference in general but a numerical integration of the equations for unequal 
initial masses with velocities of the required order has given resulting stars of 
mass-ratio 1:2. It is difficult without further investigation to say whether 
larger mass-ratios could be obtained, but this process can at any rate account for 
some differences between the masses of the components of the resulting binaries. 
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3-4. Effect of initial lateral velocities.—In the above treatment of the binary- 
forming encounter, the stars subsequently form a binary star whose components 
have highly elongated orbits. If the stars originally have slight lateral velocities 
off Ox, then the stars will follow paths which are curved but nevertheless close 
to Ox, so that the above treatment may still be expected to hold approximately 
and the resulting binary will have approximately the same energy in each case. 
However, the orbits of the resulting binaries will be elliptic or, in the extreme 
case, circular. It is important to find the maximum separation in such cases 
and this can be done by equating the energies of the different orbits. Consider 
a binary with an elongated orbit and one with a circular orbit. Let the stars be 
of equal mass m, and in the elongated orbit let the maximum separation be 2R. 
Then the energy of the binary is — Gm?/2R. In the circular orbit the energy is 


4mv* + 4mv? — 4Gm?/r, 


where visthe velocity of eachstarand r is the radius of the orbit. But v=(Gm/4r)!?, 
as can be shown from the equation of motion for either star, so the energy is 
—Gm?/4r. Equating these energies we have R=2r so that in a binary-forming 
encounter, if the stars separate to $d in the elongated case, they would separate 
to about }d in the circular case and to intermediate distances in the elliptic case. 

It is of interest to note that the angular momentum of the binary depends on 
the small initial velocities of the stars and on the inclination of these velocities 
to the line joining the stars. If inclinations up to 10° are permitted (and this seems 
reasonable), then the angular momentum can be shown to be of the same order 
as that of observed binaries. For example, if we let both stars (having the initial 
mass stated in the mth column of Tables I and II) have the maximum allowable 
initial velocities, and these are both inclined at 6° to the joining line, then the 
angular momentum of the resulting binary is 


102 x 10485 x 1015" cm? g sec"! (m= 1, 2, 3, 4), 


i.e. in the range 10° to 10°4 cm? g sec“. 

It is to be expected that the binary will interact further with the remaining 
interstellar material, with a consequent reduction in the separation of the com- 
ponents and in the eccentricity of their orbits; but it is not proposed to consider 
this subject here. 

3.5. Limitation of the elementary treatment.—The principal weakness of the 
above elementary treatment is that the mass. of accreted material is very much 
greater than the masses of the original stars. ‘This means that the gravitational 
effects of the material which is to be accreted may well be greater to start with 
than the gravitation of the stars which is alone taken into account. In fact, the 
reason for giving the elementary treatment was to show that a straightforward 
application of the accretion theory to our model does not give a reliable picture 
of the situation and to indicate in which direction improvement should be sought. 

4. A more realistic model.—We may now consider a more realistic model of 
this situation. 

We suppose that stars are condensing at various points in the star-forming 
medium which is initially at rest. We do not propose to study the reasons why a 
star should form at any one point rather than at another, nor indeed why they 
should form at all. These are problems connected with the formation of single 
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stars. We merely take for granted that stars are in fact forming around various 
“‘centres of condensation’. In a completely homogeneous medium there are 
no (local) variations in the gravitational field. It is only when condensations 
begin to occur that these local variations develop and these variations in turn 
produce movements of the interstellar material which we shall now consider. 
Let two neighbouring centres of condensation A and B be situated a distance d 
apart. The magnitude of d is to be considered as of the same order as the average 
distance between condensations in the neighbourhood, but the distance from A 
to B is to be less than the distances of either of these points from any other centre 
of condensation. We shall first briefly describe the process we envisage and then 
we shall discuss its various aspects. 

The material within the spheres EIHK and FLGJ (Fig. 1), which have centres 
at A and B, collapses under its self-gravitation to form stars at A and B. At the 
same time, the material in the cylinder EFGH collapses on to the axis AB, partly 
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under its own gravitation but mainly under the gravitational influence of the 
condensations which are formed at A and B. As the material collapses on to the 
axis, collisions between the particles of the material occur and heat is generated 
and radiated away. In this way, energy is removed from the material along the 
axis and so it remains as a column of comparatively small section extending 
between the stars at A and B. The second stage of the process then occurs in 
which the stars at A and B move towards one another under their gravitational 
attraction. In so doing they move through the column of material and sweep it 
up. We shall show in the calculations below that this column offers substantial 
resistance to the motions of the stars which can easily remove sufficient energy 
from the stars to form them into a binary. The motion of A towards B will not 
quite be rectilinear due to the disturbing gravitational influence of the neigh- 
bouring stars. This influence makes the probability very small of a collision 
between A and B. 

The first point which requires consideration is the size of the spheres ELHK 
and FLGJ. These are to be sufficiently small so that at points within one sphere 
the gravitational attraction of the material in the other is small. If the material 
inside EIHK collapsed purely under its self-gravitation, it would all arrive (as 
we shall see) at the centre A at thesametime. This will not occur in reality because 
the particles at points like I on the surface of the sphere will be retarded due to the 
attraction of the material condensing at B. The size of the sphere is arbitrary 
but it is taken so that all the material in it will arrive at A at approximately the same 
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time in spite of the retardation of some of the material due to other centres of 
condensation. The sizes of the spheres need not be the same and they will depend 
on the number, sizes and proximity of neighbouring condensations. 

A few words may be said on the behaviour of the material in the region 
EFJGHI. The process of collapse on to the axis is the same as the process 
established in the Bondi—Hoyle—Lyttleton accretion theory for the capture of 
interstellar material by stars. The difference in this case is that the column of 
material forms ahead of the stars instez. of behind as in the accretion theory. The 
main point is that if the accretion process is possible, then so is this process, since 
it uses the same assumptions about the physical condition of the medium. In 
fact the column produced in this process is more stable than that of the accretion 
theory for if a particle strays from the axis, it immediately experiences a net 
force (due to A and B) tending to bring it back. The process by which each star 
sweeps up the material in the column is precisely the mechanism of Bondi and 
Hoyle as for a single star (i.e. the star “‘ accretes ” the material). It follows from 
this that the width of the column of material lying between the stars need not 
necessarily be as small as the diameter of the stars. It also follows that the stars 
can sweep up most of the material in the column, even if (as would be the case) 
the paths of the stars do not coincide with the axis of the column. 

In considering the pair of stars A and B, we have seen how the material within 
ECFLGHK is predominantly influenced by the condensations at A and B. The 
neighbouring condensations produce a cut-off effect which determines the shape 
of the region which we have ideally represented by ECFLGHK, i.e. the region 
within which the material all collapses towards the axis AB. As regards the 
condensations and material beyond the immediate neighbourhood, we may expect 
their joint gravitational field to be negligible near AB if the condensations have 
formed from an originally uniform medium. If the medium was not originally 
uniform we may still expect this gravitational field to vary only slightly over the 
region ECFLGHK and it will therefore have a negligible effect on the processes 
outlined above. 

For the purpose of the calculations, the radii of the spheres EIHK and FLGJ 
have both been taken as }d. It is quite possible that material outside the region 
EFLGHK also is attracted and is involved in the process. It will be shown, 
however, that such material will not to any extent reduce the likelihood of binary 
formation. 

4.1. Time of collapse of a sphere of material.—Our first calculation is to deter- 
mine the time the material in the sphere EIHK takes to collapse assuming a free 
fall of the material under its self-gravitation. If we consider a particle originally 
at rest at a distance R from A, then it is attracted towards A by a force which is 
the same as if all the material within a sphere of radius R were concentrated at A. 
This force, when the particle is at a distance ry from A, is $GrR%p/r?. By 
integrating the equation of motion we obtain the velocity v of the particle, 

hv? = $GrR%o(R/r —1), 
then the total time for the fall of the particle to A is 
‘R 
| =(3n/32%Gp)2 
0 
=0°543(Gp) **, (11) 
which is independent of R, so that all the particles reach A at the same time. 
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In our initial statement of the process it was implied that the stars condense 
at A and B and then begin to attract each other. In fact, of course, the attraction 
begins as soon as any sort of condensation occurs at A and B. We can show, 
however, that the condensations move a negligible distance in the time (11). 
The distance moved from rest by A in time t is } ft®, where f is the acceleration. If 
all the material in the sphere FLGJ is considered concentrated at B, then 
f = $Gn(d/4)%p/d?, and putting t equal to (11), we have }/t? = 0-0096d. 

4.2. Time of collapse of cylindrical region.—In order to obtain an estimate of 
the time required for the region EFJGHI to collapse we shall calculate the time 
required for a particle at C to fall to the centre O under the influence of the 
material in the spheres EIHK and FLGJ on the assumption that this material 
attracts as if it were all concentrated at A and B respectively. At any time ¢ let y 
be the distance from AB of the particle, originally at C, and let v be its velocity. 
The acceleration of the particle at this time, due to the attraction of A and B, is 


2G§n(td iy 
ed? OU 


Integration of the equation of motion gives 
Jo? = 2G. §a(bd)Ppl (4d? + y*) 1 (ha? + yo) "4 


Yo being the initial value of y, i.e. }d. The time for the particle to fall from 
C to Ois 





th dy 


o Uv 
This was evaluated numerically and was found to be 


1-74(Gp)*?. (12) 
This must be compared with the time required for the star A or B to fall to O. 
We can obtain a lower limit of this by assuming that these stars experience no 
resistance. It can easily be shown that the time required for two particles each 
of mass $7(}d)%p to come together from rest at a separation d is 


(377)"(Gp)-"? = 3:07(Go)-™”. (13) 
So, as (12) is little more than half (13), it is seen that the cylindrical region 
EFJGHI has adequate time to collapse before the separation of the stars is 
appreciably reduced. If we had taken the self-gravitation of the cylindrical 
region into account, the time of collapse would have been even shorter. 

4.3. The approach of the stars.—We must now consider what happens when 
the stars approach one another and encounter the column of material. We shall 
assume that the stars move under their mutual attraction and that they capture 
all the material in the column as they move through it. To start with, we shall 
assume that the material is uniformly distributed in the column which extends 
from A to B. We shall suppose the mass of material per unit length of column 
isa. Then considering the axis Ox to be coincident with AB, the coordinate of B 
is +}d. Writing down the equations for the star B we obtain equation (2) as 
in Section 2. We also have 


dm) | Gnt dm, 
ge ieee" 
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Since v is negative 
dm 
dx 


and so m=—ax+M, (14) 


where M is the mass of the star when it reaches the origin. ‘The equation for 
v® becomes 
G(M — ax) 
is 

of which the solution is 

v= Hl gs + 3M?alnx — 3Ma*x + jars?) +4 | (15) 

2(M—ax)L\ x , 

the constant A being determined by the fact that v=o when x=}d. Now the 
total energy (kinetic + potential) of the two stars at any instant is 


(M — ax)v? — G(M — ax)?/2x, 
and using (15) this becomes 


G 

2(M—ax) 
In this, as xo, the logarithmic term predominates and becomes large and 
negative. Consequently, the total energy of the pair of stars becomes greatly 
reduced so that they cannot again separate. If the stars were to continue to 
collect material until they get right to O, then In x — co and the stars would be 
unable to separate at all. Furthermore, this would occur however small is a. 
We suppose in reality that the process breaks down before O is reached. This could 
be caused by the stars being slightly deflected off Ox by the gravitational field of 
neighbouring stars. To take a particular example, suppose the process discon- 
tinues when the stars are 4d apart, then their subsequent separation would not 
exceed 039d. (In making this calculation, we used a= M/d, which is what it 
would be if the material in EFJGHI were uniformly distributed along AB.) 

The above analysis is based on the assumption that the material in the column 
remains at rest until it is collected by one of the stars. This assumption is invalid 
for two reasons. In the first place, when the material first collapses on to the 
axis, most of it (e.g. particle D in Fig. 1) will have a velocity component along 
the axis away from O. The effect of this is that the stars will be encountering 
material moving against them. They will therefore lose more momentum than 
in our treatment and consequently this treatment does not overestimate the 
resistive effect of the material column. In the second place, we have ignored the 
gravitational interaction between the stars and the unaccreted parts of the column. 
This is a difficult factor to assess, but it seems to be an effect which would largely 
cancel itself. For, consider the star A: at the beginning of its approach towards 
O, it will be attracted by the material in the column and this will cause its 
momentum to be increased. However, the material in the column (to the left 
of QO) will be attracted by A and will begin to move towards A. Consequently, 
when the material collides with A, the momentum which A originally gained 
will largely be lost in destroying the momentum of the material. 

A close examination of the above analysis shows that the occurrence of the 
logarithmic term depends on the assumption of a uniform column. However, 


[3.M2a In x + $a°x? — 6Ma*x + 3M?a+ A]. 
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this is not essential to the process; all that is necessary is that material is captured 
by the stars sufficiently close to the origin O. This explains why material from 
outside EFLGHK cannot appreciably affect the process. The physical principle 
involved is simple: it is that if a mass L moving with velocity U collides with and 
captures a mass / at rest, then the energy loss is }L/U?/(L+1/), which is large 
if Uis large. The velocities of the stars oo as they approach O. 

4.4. Comparison with the previous process.—In comparison with the binary- 
forming encounter considered in Sections 2 and 3, in this process there is no limit 
on the initial masses of the stars which condense at A and B. Also, the mass of 
material subsequently captured is of the same order as the initial masses of the stars 
so that the possibility of a star having a second chance to form a binary is not 
ruled out. Nevertheless it seems likely that a process of the type here con- 
sidered could only have been important at a very early stage in the evolution 
of the stars. In this process it seems quite possible that any reasonable mass 
ratio of the binary components could be obtained. Although more of the 
interstellar material is used up in the process, a good percentage is left for 
further effects. 

5. Conclusion.—In putting forward this mechanism for binary star formation 
it is hoped that it will indicate a way by which binary stars can be formed from the 
‘star-forming medium without the postulation of very close initial condensations. 


It remains for me to express my thanks to Professor W. H. McCrea for his 
help and encouragement during the development of the work contained in this 
paper. The idea of binary stars being produced as a result of the resisted motion 
of single stars through interstellar material was first conceived by Professor 
McCrea and it was this idea that led to the present work. 


Royal Holloway College, 
Englefield Green, Surrey : 
1954 June 28. 
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TOPOCENTRIC LIBRATIONS FOR AN ECLIPSE STATION 


C. A. Murray 


(Communicated by the Astronomer Royal) 


(Received 1954 November 12) 


Summary 


The selenographic coordinates of an observer who is instantaneously on 
the axis of the Moon’s shadow during a solar eclipse are diametrically 
opposite to those of the Sun. The latter, and hence the former, may be 
obtained directly from the Nautical Almanac. The present paper gives 
a short method for computing differential corrections to these coordinates, 
giving the topocentric librations for any eclipse station. "The method makes 
use of computations which will have been made in any case for the prediction 
of the local circumstances of an eclipse. 





In a recent paper (1), Atkinson has given a simple method for computing the- 
topocentric librations and the position angle of the Moon’s axis appropriate to 
any lunar observation. His method is to apply parallactic corrections to the 
geocentric librations (‘‘ Earth’s selenographic longitude and latitude’) which 
are tabulated in the Nautical Almanac. 

In the case of solar eclipse observations the computation may be shortened by 
making use of the selenographic coordinates of the Sun instead of those of the 
Earth. The former are readily obtainable from the Almanac by linear interpola- 
tion, whereas the latter require second differences at least. For an observer who 
is instantaneously on the axis of the shadow, i.e. one for whom the eclipse is 
central at the instant of the observation, the topocentric librations are dia- 
metrically opposite to the selenographic longitude and latitude of the Sun. 
Accordingly if /,, 6) are the topocentric librations for this observer, 

1, = 270° — Sun’s selenographic colongitude, 
b= — Sun’s selenographic latitude. 

‘The topocentric position angle of the Moon’s axis may be obtained from the 
general formulae, given in the N.A. “ Explanation ” from 1931 to 1941, modified 
to take account of the circumstances peculiar to eclipses. The most suitable 
form of the general formula for C’, the position angle of the Moon’s axis without 
physical libration, is 

sin C’ = —sin 7 cos (a — §%’) sec b’, (1) 
where « is the topocentric R.A. of the Moon, b’ the optical libration in latitude, 
and 7, §%’, are mean elements of the Moon’s equator tabulated in the Almanac at 
ten-day intervals. No sensible error is caused by writing sec b’=1 in computing 
Cy’, the value of C’ for the observer on the axis of the shadow. Also, for this 
observer, 

a = Greenwich sidereal time — p, 
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‘where p, the Greenwich hour angle of the point on the celestial sphere toward 
which the axis of the shadow is directed, is tabulated at ten-minute intervals 
throughout the eclipse. Accordingly we compute C,’ from 
sin Cy’ = —sin i cos (G.S.T. —p— 8’). (2) 
The formula for the physical libration of the Moon’s axis reduces in the 
present case to 
C” = —0°-020 cos (I” — R+ 21’) F 0°-003, (3) 
since the libration in latitude must be nearly zero and the Moon near one of its 
nodes. I’, are mean elements of the Moon’s orbit tabulated at ten-day intervals 
in the Almanac and I’ is the optical libration in longitude. It is sufficient to sub- 
stitute /, for /’ in (3). The upper or lower sign in (3) is to be taken according as 
the Moon is near the ascending or descending node. The topocentric position 
angle of the Moon’s axis for the observer on the axis of the shadow is then 
Cy=C,'+C", (4) 
where C,’ and C” are computed from (2) and (3) respectively. 
The topocentric librations and position angle of the Moon’s axis for an 


observer who is situated at any point where the eclipse is visible may be obtained 
by the application of differential corrections to J), by and C5. 


Let l=)+Al, 
b=b,+Ad’ 
be the topocentric librations for this general observer; (A/’, Ab’) are then the 
coordinates of the centre of the apparent lunar disk, with respect to the point 
on the lunar surface which is on the axis of the shadow (i.e. the centre of the 
apparent disk for the observer on the axis of the shadow). 
(z—¢) 


B 
O Ww +: 








OC’ 
Fic. 1. 

In Fig. 1, let O be the observer on the axis of the shadow, O’ the general 
observer, A the Moon’s centre, and B the intersection of the axis of the shadow 
OA with the plane through O’ perpendicular to OA. In the usual Besselian 
coordinate system and notation, 

x—&=m sin M, 

y-—n=m cos M (5) 
are the coordinates of B with respect to O’, and the distance AB is (z—). Let 
w be the angle O’AO, then 


tan w= —y (6) 





and the total change in optical libration in moving from O to O’ is w in position 
angle /+180°. With an error not exceeding 4 parts in 10* we may write 


3 = cosec 7 
( 
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where 7, is strictly the Moon’s equatorial horizontal parallax at the instant of 
observation, but may be taken as the value for the time of conjunction in R.A. 
of the Sun and Moon which is given under “ Elements of the Eclipse” in the 
Almanac. Also, since z~60, |{|<1 and m <o-6 we may write (6) as 


w= mr ¢ (7) 
with a maximum error of 0°-o1. This is sufficiently accurate for all purposes at. 
the present time. (See Atkinson, Joc. cit. p. 452.) 

Since the libration in latitude is nearly zero the extremities of the Moon’s 
axis of rotation are for all practical purposes on the apparent limb. 


N 








Fic. 2. 


In Fig. 2, which represents the topocentric view of the Moon for the observer 
O, P is the north pole of the Moon, W the west limb, N the direction of the north 
celestial pole, and D, D’ the centres of the apparent disks for the observers O, O’ 
respectively. Then evidently, 


Al'= +w sin (M—C,), (8) 
Ab’ = —w cos (M—C,), 
since selenographic longitudes are measured positively towards the west limb 
and latitudes towards the north. Combining (5) and (7) with (8) we have 


Al’ = +7 {(x—£) cos Co—(y—y) sin Co}, 
Ab’ = —m_ {(x—£) sin Cy+(y—) cos Co}. 


The change in the topocentric R.A. of the Moon in moving from O to O’ is 
evidently 


(9) 


Aa =m, (x—€) sec d, 


where d is the declination of the point on the celestial sphere toward which the 
axis of the shadow is directed. Hence, logarithmic differentiation of (1) gives 
the change in the topocentric position angle of the Moon’s axis 


AC’ = —7¢(%—£)sec d tan(G.S.T.—p— 9’) tan C;,’ (10) 
neglecting a small term in Ad’ tan b’. 
The main disadvantage of the above method is the necessity for computing 


C, from first principles. If this were tabulated in the A/manac at intervals during 
the eclipse the work would be considerably shortened. 
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The method has advantages over Atkinson’s method even when C, has to be 
computed. It makes use of quantities derived in the course of the prediction of 
the local circumstances, instead of requiring the computation of the apparent 
zenith distance, parallax and parallactic angle of the Moon. Although it requires 
the same number of interpolations from the Nautical Almanac (in addition 
to those needed in the course of the prediction), these are all linear, whereas at 
least three of the quantities in Atkinson’s method, the geocentric librations and 
position angle of the Moon’s axis, require second differences. 


In conclusion I wish to thank Mr D. H. Sadler for his helpful criticisms of 
this paper. 
Royal Greenwich Observatory, 
Herstmonceux Castle, 
Hailsham, 
Sussex : 
1954 November ro. 
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AROUND h AND x PERSEI 
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(Received 1954 November 29) 


Summary 


Photoelectric colours and magnitudes for thirty early-type stars, classified 
as supergiants belonging to the double cluster in Perseus, have been determined 
with the aid of the Lick 12-inch refractor. Measures of four stars in Harvard 
Standard Region C12 were used for the tie-in to the (P, V) system. 

A general increase of absorption is found as one passes from h to y Persei, 
in accordance with the earlier conclusions of Miss Pismis. A considerable 
non-uniformity of interstellar absorption exists around the cluster. 

The supergiants of type B2 and the Bo, Bo-s5, B1 and B2 stars of 
luminosity class III yield a distance modulus of 11™-7 for the cluster of 
supergiants, corresponding to a distance of 2-2 kps. Using the new value of 
the distance modulus, absolute magnitudes have been calculated for super- 
giants of other spectral classes. The number of stars used in this investigation 
is too few to give an estimate of the dispersion in M, for these early-type 
supergiants. 





Highly luminous early-type stars are generally too far away from the Sun to 
give a reliable determination of their absolute magnitudes from parallax data. 
We therefore resort to the study of magnitudes and colours of objects of this kind 
that form part of large stellar aggregates with known distance moduli. The 
double cluster in Perseus, the Orion aggregate and NGC 6231 have, more or less, 
a fair sample of such supergiants. The double cluster in Perseus, however, con- 
tains a larger number of B and A type supergiants than the other two. It also 
contains a few very luminous stars of type M. It is natural, therefore, that our 
knowledge of the absolute magnitudes of these objects, and their dispersion in 
luminosity, should rest heavily on the results that can be derived from this cluster. 

Much of our knowledge concerning the remarkable cluster of supergiants 
around h and yx Persei depends on a thorough investigation by Bidelman (1). 
He employed the MKK system of classification and separated the likely cluster 
members from the field supergiants. More recently Morgan, Whitford and 
Code (2) have provided more accurate classes on the MK system for some of the 
cluster members sorted out by Bidelman. 

It is the purpose of this investigation to determine accurate magnitudes and 
colours for those supergiants belonging to h and y Persei which can be “ resolved ” 
by a photoelectric photometer on a visual refractor of medium focal length. 

The photoelectric observations of colours and magnitudes for 30 supergiants 
in and around the double cluster were carried out on four nights in 1952 September 
using the 12-inch refractor of Lick Observatory. The photometer employed a 
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1P21 photomultiplier tube, the output of which was amplified by a d.c. amplifier, 
designed by Kron, and fed into a Brown recorder. The filters used were Corning 
5562 and 3385 giving effective wave-lengths of approximately 4200 A and 5200A 
respectively. ‘The deflections were obtained in the order blue, yellow, blue, 
followed by readings on the sky in the two colours. Also, deflections caused by a 
standard source provided a check on the stability of the overall response of 
‘the system. Mean extinction coefficients determined from Mount Hamilton 
extinction data gathered over several years were utilized in the reductions. 
Nightly extinction coefficients were also derived by Eggen’s method (3) for three 
of the four nights, and these proved to be very close to the mean coefficients used. 
The reductions were carried out using the mean extinction coefficients as well as 
the values determined observationally for the night under consideration. 
No sensibly different values were obtained for magnitudes or colours, and hence 
the mean extinction values were used throughout finally. 

The instrumental magnitudes m, and colours C;,, were transformed to the 
(P, V) system by means of nightly comparisons with stars A, C, D, F of Eggen’s 
photoelectric sequence (4) in the Harvard Standard Region C12. While the V 
magnitudes used were those of Eggen, the P— V values used were revised values, 
kindly supplied by Dr Kron in advance of publication. The transfers were made 
every night when C12 and the double cluster had nearly the same zenith distance. 
The transformation equations used are the following :— 


P—V=A+BC,, (1) 
V=m,+D+F(P-V), (2) 
A=0°417+0°012, 

B=1-055 +0°005, 
D = 6'494 + 0°003, 
F=—0-116+0:006. 


From the individual values of V and P—V obtained each night for the stars 
under observation the probable errors for one observation were found to be 


P-—V : p.e.= +0™-009, 
V: p.e.= +0™-012. 


In Table I are listed the relevant data for thestars included in this investigation. 
The first column contains the HD number, the BD number being given for one 
star not listed in the Henry Draper Catalogue. The second and third columns 
list the observed magnitudes V and colours P— V, these being means of three or 
four observations (rarely only two) on separate nights. Column 4 contains 
information on the spectral types; those with asterisks are due to Bidelman, 
the rest have been recently classified by Morgan (2) on the MK system. Finally 
in column 5 we have the colour excesses derived from the intrinsic colours of 
unreddened stars, on the (P, V) scale, as transformed from the (B, V) values given 
by Morgan, Harris and Johnson (5). 

The available sources of reliable magnitudes or colours in the h and x Persei 
region are the investigations of Wallenquist (6), Oosterhoff (7), and Stebbins, 
Huffer and Whitford (8). Fourteen stars in Table I have also been included in 
Wallenquist’s photoelectric studies. Assuming that Wallenquist’s system of 
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TABLE I 
Colours and magnitudes of early-type supergiants around the double cluster 





HD Number or 


BD Number V P—V _ Spectraltype Colour excess Remarks 





12856 +0°08 *Bo-5 ne 0°56 MWC 25 

12953 +0°45 Arla 0°59 

13051 8-75 —o'ol *Bo-s 0°47 MWC 27 

13267 6°42 +0°17 Bs Ia 0°52 

13476 6°49 +0°46 A3 Iab 0°54 

13744 7°65 +0°60 Ao Iab 0°77 

13621 8°17 —0'09 *Bo-s 0°57 ADS 1714 

13745 7°94 +o-o1 Bo III o's 

13841 7°46 +0°08 B2 Ib 0°52 

13854 6°54 +0'13 Br Iab 0°59 MWC 31 
56° 473 9°13 +0°09 *Bin 0°55 MWC 441 

13866 7°57 +0°05 B2 Ib 0°49 

13900 9°26 +002 *Bin 0°48 

13969 8-86 +o°15 *Bi o'61 

14250 g'01 +0'20 *Ban 

14322 6°82 +018 B8 Ib 

14321 9°31 +015 *B2 

14422 8-96 +0°36 *Bine ° MWC 37 

14433 6°45 +0°42 Ar Ia 

14434 8°57 +0°04 *O7n 

14476 8°79 +0°25 *Bo-s 

14520 9°23 +0°10 *Br ° ADS 1810 

14535 7°56 +0°57 Ar la 

14542 7°03 +0°47 B8 Ia 

14818 6-28 +o°15 B2 Ia 

14899 7°47 +-0°29 *Ao 

14956 7°26 +0°58 B2 Ia 

15316 7°29 +064 A3 Iab 

15497 7°08 +0°63 B6 Ia 0°96 

14143 6-70 +0°34 B2 Ia 0°78 





* Spectral classification by Bidelman. Those not starred are due to Morgan (2). 


magnitudes and colours can be converted to the (P, V) system by a linear trans- 
formation, we get the following relations, by least squares, using stars fainter than 
g™-o for the purpose: 
P=m,, + 0°037 —0°140 (P—V), (3) 
P—V=o'119+0'966 C,,. (4) 
A similar procedure for Oosterhoff’s photographic magnitudes gives the relation 
P=m,, +0051 +0:163 (P—V). (5) 
The values derived from equations (3), (4) and (5) are compared with those from 
Table I in Fig. 1. The central diagram of Fig. 1 shows that equation (4) is a 
good approximation, as the colour differences between the two photoelectric 
systems are negligibly small over the range of magnitudes observed. Below 
magnitude g-o a randomness sets in in the measures that might be considered 
to indicate the rapid approach of the limiting magnitude for accurate measurement 
with Wallenquist’s experimental arrangement. The two magnitude systems, on 
the other hand, are not completely described by equations (3) and (5). The 
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similarity between the trends in the upper and lower diagrams of Fig. 1 is probably 
due to the heavy weighting of Wallenquist’s scale and zero point with the aid of 
Oosterhoff’s magnitudes. 

Twenty-two stars in Table I have colours on the C, scale measured by 
Stebbins, Huffer and Whitford. These colours can be transformed to the (P, V) 
system by the linear relation 


P—V=0-199+ 1-722 Cj. (6) 


There is, however, a large scatter among the points. A better fit is obtained 
by the curve 


P—V=0-208 + 1°833 C,— 1-042 C,?. (7) 





























Fic. 1.—Comparison of the photoelectric observations with the results of Wallenquist and Oosterhoff. 
Pw, (P—V)w and Po are calculated values using equations (3), (4) and (5) respectively. 


A plot of observed colours against spectral type, as seen in Fig. 2, shows the 
non-uniformity of interstellar absorption around the double cluster. Fig. 3 
shows the variation of colour excess (100 colour excess is plotted) over the region 
containing the two clusters. There is a general increase of absorption as one 
passes from h to x Persei, in agreement with the results of Miss Pismis (9); and 
superimposed on this is the irregular variation caused by a patchy distribution of 
the absorbing clouds. Consequently, any observed colour-magnitude array 
would show a large dispersion in values along any of the sequences of interest. 

In Fig. 4, I have plotted apparent visual magnitudes corrected for interstellar 
absorption (2-8x colour excess) against spectral type for the stars listed in Table I. 
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Fic. 2.—Colour-spectrum diagram for supergiants around the double cluster. 














° 


l 


Fic. 3.—Variation of colour excesses in the region around h and x Persei. 
The figures shown are 100 X colour excess. 
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The solid dots represent stars classified by Morgan, while the circles indicate 
the remaining stars having Bidelman spectral types. The trend followed by the 
stars of luminosity class I clearly shows that it may be reasonable to assume some 
of the stars classified by Bidelman to fall in the category of luminosity class III. 
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Fic. 4—V magnitudes, corrected for interstellar absorption, for supergiants around hand x Persei. 
Filled circles are for stars classified by Morgan, while the open circles represent stars having 
Bidelman spectral classes. 


Adopting this assumption, and using the absolute magnitudes given by Keenan 
and Morgan (10) corrected by —o™-1 to bring them in accordance with the visual 
magnitudes used here, I get the values for the distance moduli given in Table IT. 
The distance moduli derived from the B2 Ia and Bz Ib stars agree remarkably 
well with those of Bo, Bo-5, B1 and Bz stars of luminosity class III. Assuming 
then that the distance modulus of the B-star aggregate is 11™-7, we get the values 
of absolute magnitudes given in Table III for the other spectral types. The 
quantity in parentheses after each value of M, indicates the number of stars used 
for its determination. It becomes apparent from Table III as well as from Fig. 4 
that the absolute magnitudes of the super-giants of class Ia do not possess the 


Taste II 
Determination of distance modulus of the cluster of supergiants around h and y Persei 





Spectral type My M, m,—M 





Bo-s III 6°95 —4°5 11°45 
Br Ill 7°40 —4°4 11°80 
B2 Ill 7°44 —4'2 11°64 
B2 Ia 4°51 —7'! 11°61 
B2 Ib 6-10 —5'8 I1°gO 
Mean m,— M=11°7 
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same values for all spectral types. A decrease in M, from Bz to Az is quite 
noticeable. This cluster contains very few stars to permit a reliable study of the 
dispersion in M, among the early-type supergiants. Nevertheless it does appear 
that the intrinsic dispersion may be much smaller than is generally assumed. 
If the cluster of supergiants is spherical and has an angular diameter of four 
degrees, the depth of the cluster in the line of sight will cause a scatter of o™-15 
in the values of M,. This value is on the assumption that the intrinsic dispersion 
experienced by the supergiants is zero. 








TasBLe III 
Derived absolute magnitudes for early-type supergiants 
Spectral type Ib Ia 
Bs —6°7 (1) 
B6 —7°3 (1) 
B8 —6:1 (1) —6°7 (1) 
Al —6°5 (3) 





A very accurate luminosity classification of many more supergiants around 
hand x Persei, along with photoelectric magnitudes and colours, will be required 
to determine the dispersion in M,. The ideal regions for getting more information 
about absolute magnitudes of all representative supergiants are the Magellanic 
Clouds. It is natural to hope, therefore, that those having access to these external 
galaxies would be able to solve the problem before long. 


I am grateful to Dr C. D. Shane and Dr Akbar Ali for having generously pro- 
vided all facilities for work at the Lick and Nizamiah Observatories respectively. 
It is a pleasure to acknowledge indebtedness to Drs Kron, Bok and Miczaika for 
critical discussions, and to Mr M. Ramanathan of the Nizamiah Observatory for 
having carried out a major portion of the reductions. 


Nizamiah Observatory, 
Begampet, Hyderabad, 
Deccan, India : 
1954 October 1. 
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MAGNITUDES AND COLOURS OF SOME MEMBERS OF THE 
PERSEUS CLUSTER 


M. K. Vainu Bappu 
(Communicated by the Director, Nizamiah Observatory) 


(Received 1954 November 29) 


Summary 


Photoelectric colours and magnitudes have been determined for twenty- 
five stars considered to be members of the moving cluster in Perseus. 
A distance modulus of 6™-83, corresponding to a distance of 230 ps, has 
been derived. This value is nearly twice the value of 120 ps derived by Miss 
Roman and Morgan from the upsilon component of proper motion. 





In a recent investigation, Miss Roman and Morgan (1) have examined the 
visual clustering around the F5 supergiant « Persei. They have shown that 
twenty-five stars of types Bz to A2 appear to share the space motion of the super- 
giant, and that a fair correlation exists between spectral type and apparent 
magnitude. The evidence provided by Roman and Morgan seems to favour the 
reality of this aggregate being a moving cluster. 

Earlier, Smart and Ali (2) examined a list of 48 stars spread out over a large 
area of the sky and which had been ascribed tothe moving cluster. They concluded 
“ (a) that a considerable proportion of the stars in the so-called Perseus cluster 
cannot belong to a moving cluster, and (4) that a nucleus of stars remains which 
may possibly be regarded as a cluster formation”. Roman and Morgan comment 
that if « Persei, or some B stars of the visual clustering around the supergiant, 
had been considered as representatives of the group by Smart and Ali, the reality 
of this smaller group would have been easily noticed. 

To settle the question of the existence of the Perseus moving cluster definitely, 
one weuld need reliable radial velocities and a homogeneous set of proper motions, 
which do not exist at present. Establishment of an accurate colour-magnitude 
array would decide whether the stars form part of a cluster. Whether they are 
members of a moving cluster can be determined only if a distance modulus 
(easily derivable from the colour-magnitude array) agrees, within the usual limits 
of tolerance, with that determined dynamically. In this paper a magnitude- 
spectral type diagram plotted from photoelectric data gives us a distance modulus 
for the cluster. 

Twenty-five stars that were shown to be definite cluster members by Roman 
and Morgan have been included in this study. The photoelectric observations 
were made on three nights in 1952 September with the 12-inch refractor of Lick 
Observatory. The observational procedure and method of reduction are similar 
to those described in an earlier paper (3). Thirteen stars of this list have only one 
observation each. The consistency of the P—V and V values derived for the 
remaining twelve stars indicates that there is no lack of observational accuracy in 
the colour and magnitude values derived from a single observation. One source 
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of error, however, may exist, and that is the recording of the sensitivity of the 
amplifier for the single observation. The Lick 12-inch photometer has gain- 
switch positions which change in steps of half a magnitude and, consequently, a 
mistake in recording will tend to alter the star’s V value by +0™-5, while P—V 
value remains unchanged. The magnitudes obtained in this paper agree with 
those derived by Roman and Morgan without any serious discordances; thus 
such a source of error can scarcely exist. The tie-in with the (P, V) system was 
effected by nightly comparisons with C12, where the V values were those given 
by Eggen (4) and the P—V values were revised values kindly supplied by 
Dr Kron in advance of publication. ‘line probable errors for a final value of V 
and P—YV are the following: 
P—V : p.e.= +0™-008, 
V : p.e.= +0™011. 


Table I gives the magnitudes and colours of the cluster members investigated. 
The colour excesses have been computed by using the intrinsic B—V colours of 


TABLE I 
Magnitudes and colours of members of the moving cluster in Perseus 
. 7 No. of Spectral 
Star V P-—V che. | Ge E | oe M 
HD 20365 5°20 —o'20 2 B3 V o"19 4°69 —2°'1 
HD 20418 5°08 —o'19 2 Bs V 0°16 4°63 —2°2 
HD 20391 8-o1 0700 3 A1V O14 7°62 +o°8 
HR ro11 5°38 —o'20 2 Bs V O's 4°96 —1°8 
HR 1029 6°13 —o'2I 2 B6 V o'12 5°79 —1'0 
HR 1034 5°03 —o'22 3 B3 V O'17 4°55 —2°'2 
HD 21375 7°55 —0'03 3 A1V orl! 7°24 +04 
HD 21479 7°36 —o'05 2 A2V 0°06 7°19 +04 
34 Per 4°72 —0o'24 2 B3 IV O's 4°30 —2°5 
HR 1037 5°65 —o19 2 B6 V O14 5°26 —1.5 
HD 21481 7°70 —o'00 «2 Ao V O17 7°22 +0°4 
HR 1051 5°87 —o17. 2 B8 IV o*10 5°59 —1°'2 
HD 20701 8-43 —0'03 I A1V o'll 8-12 +1°3 
HD 22401 7°50 —o12 I Ao V 0°05 7°36 +06 
HD 20961 7°74 —0'03 I AoV O°14 7°35 +06 
HD 21091 7°60 —orll I Ao V 0°06 7°43 +06 
HD 21181 6-88 —o'l3 I Bg V Ovll 6°57 —o'2 
HD 21398 7°44 —o'l2 I Bo V o'12 7°10 +03 
HR 1047 6°27 —o'ol I Bs V 0°34 5°32 —I1°5 
HD 21641 6°77 —0'l3 I Bg V o'll 6°46 —0°3 
HD 21672 6°61 —o'16 I B8 V orll 6°30 —0'5 
HR 1063 5°53 —o'22 I B8 III 0°05 5°39 —1"4 
HD 21931 7°41 —or'll I Bg V 0°13 7°05 +03 
x Per 4°27 —0'20 I Bse Orls 3°85 —2'9 
HD 22136 6-92 —o'ls I B8 V o°12 6°58 —o'2 


Morgan, Harris and Johnson (5), converted to the (P—V) scale. Fig. 1 is a plot 
of spectral type against observed colour. The absorption is non-uniform. In 
particular HR 1047 is highly reddened. But its membership in the cluster can 
scarcely be disputed, as will be seen from the colour-magnitude array in Fig. 2. 
Fig. 3 is a plot of spectral type against V, now corrected for absorption by sub- 
tracting 2-8 x the colour excess (see col. 7 of Table I). The scatter among the 
stars in Fig. 3 is what we would expect for the main-sequence stars of early type. 
Most of the stars, therefore, do form a cluster. 
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Fic. 1.—Colour-spectral type diagram for members of the Perseus cluster. Open circles are intrinsic 
colours, while the observed colours are represented by filled circles. 
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Fic. 2.—Colour-magnitude array for members of the Perseus cluster. 
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The distance moduli of the stars were next derived from the V magnitudes, 
corrected for absorption, and the absolute magnitudes of Keenan and Morgan (6). 
These absolute magnitudes have been corrected by —o™-1 to bring them into 
accordance with the visual magnitudes used here. The mean distance modulus 
obtained from twenty stars of luminosity class V, with the exception of HD 20701, 
is 6™-83 (239 ps). HR 1063, which is classified as B8 III, has a distance modulus 
of 8™-4. If the classification of the star as a giant is correct, it cannot be con- 
sidered a cluster member. The last column of Table I gives individual absolute 
magnitudes for the stars using the distance modulus 6-83. 
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Fic. 3.—Corrected visual magnitude-spectral type diagram for members of the Perseus cluster. 





The distance modulus obtained for the cluster by Roman and Morgan from 
dynamical considerations is 5™-4 (120 ps). The distance derived photometrically 
by me is thus nearly twice that derived from the proper motion data. It may be 
noted that Roman and Morgan derive the distance of 120 ps from the upsilon 
component of proper motion. The proximity of the convergent point to the solar 
antapex might very well introduce a large error in any distance derived by this 
means. 

An interesting anomaly is the parallax of « Persei. The trigonometric parallax 
is 0”-029 +0”-005 (7) and using the corrected visual magnitude of Roman and 
Morgan, we get an absolute magnitude of — 1-3, a value incompatible with its 
spectral characteristics. On the other hand, if we use the distance modulus as 
obtained in this study, the apparent magnitude, corrected for absorption, would 
have the low value of magnitude 2:2. 
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I am indebted to the Directors of the Lick and Nizamiah Observatories for 
providing me with all necessary facilities at their respective institutions, and to 
Mr M. Ramanathan for his help in the reduction of the observations. 


Nizamiah Observatory, 
Begampet, Hyderabad, 
Deccan, India : 
1954 October 5. 
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ENERGY LEVELS IN Caxv 


David Layzer 
(Communicated by Donald H. Menzel) 


(Received 1954 December 3) 


Summary 


A discussion of the structure of term spectra in isoelectronic sequences. 
shows that spectral terms fall naturally into certain groups, which we call 
** complexes’. A complex is completely specified by the set of principal 
quantum numbers (m) and the parity p, and contains one or more entire 
configurations. The interval between terms that belong to the same complex 
varies asymptotically as the nuclear charge Z for large values of Z, while 
the interval between terms that belong to complexes with different sets of 
principal quantum numbers () varies asymptotically as Z*. The electrostatic 
coupling between configurations whose terms belong to different complexes 
tends to vanish with increasing Z; but the coupling between configurations. 
whose terms belong to the same complex approaches constancy. This 
circumstance greatly simplifies the discussion of configuration mixing in 
highly ionized atoms. 

In Caxv the configurations (1s) (2s?) (2p?) and (1s*) (2p*) are strongly 
coupled, but the complex (122‘) is weakly coupled with higher complexes. 
On the approximation of zero coupling between (122*) and higher complexes, 
the energy levels of (172*) are the eigenvalues of a tenth-order matrix, which 
reduces to two fourth-order matrices and a second-order matrix. We have 
evaluated the matrix elements partly by extrapolation along the isoelectronic 
sequence and partly by theoretical considerations. The predicted wave- 
lengths for the transitions 2p? °P,_,,, 9P;..9 are 5456 A and 5650 A respectively. 
These agree with the observed wave-lengths of the coronal lines AA 5445, 56904 
to within the uncertainty of the calculation, thus confirming the identifications 
of Edlén and Waldmeier. 





1. Introduction 

In 1942 Edlén (1) published an epoch-making paper in which he identified 
twenty-three coronal emission lines with transitions from metastable levels in the 
ground configurations of very highly ionized atoms. Only two of the identifica- 
tions were listed as doubtful. One of these, the identification of the yellow line 
A 5694 with the magnetic-dipole transition p?°P,,, in Ca xv, has aroused con- 
siderable controversy and stimulated several investigations because of the fact 
that no ion definitely known to produce coronal emission lines has an ionization 
energy approaching that of Ca xv(814eV). 

In 1947 Waldmeier (2) discovered a faint emission line in the coronal spectrum 
at A 5445, which he subsequently (3) identified with the magnetic-dipole transition 
p?*P,.,, in Ca xv. 

Shklovsky (4) rejected the identifications of AA 5694, 5445 on the ground that 
the intensity ratio /(5694)//(5445)=6 observed by Waldmeier is incompatible 
with any known mechanism for exciting the levels 2p?*P,, , in Ca xv. However, 
Mme Pecker, Billings and Roberts (§) have since discovered that A 5445 is blended 
with a Fraunhofer line. Making allowance for the blending, they derived an 
intensity ratio of 1-6, in satisfactory agreement with theoretical expectations. 
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These writers also measured the Doppler widths of the red (Fe x, 46374), 
green (Fe xiv, 4 5303), and yellow (A 5694) lines on two occasions when all three 
lines appeared simultaneously in the same (projected) region of the corona. By 
assuming that the three lines originated on each occasion in a single region, 
characterized by a unique kinetic temperature, they derived the atomic weight 
40+2 for the ion producing the yellow line. This result tends to confirm the 
identifications of Edlén and Waldmeier, though it does not exclude the possibility 
that argon (atomic weight 39-9) may be responsible for the yellow line. 

The most serious objection to the identifications of Edlén and Waldmeier 
rests on a very thorough semi-empirical study by Garstang (6) of the configuration 
(15?) (2s?) (2p?) in the isoelectronic sequence headed by C1. On the theory em- 
ployed by Garstang, the five energy levels of this configuration are expressed in 
terms of five atomic parameters. Garstang determined the values of these 
parameters empirically in the first ten members of the isoelectronic sequence, for 
which the experimental data are nearly complete. He then (7) extrapolated these 
values to Ca xv, and used the extrapolated parameters to calculate the positions 
of the energy levels. The wave-lengths predicted for the transitions °P,_,, 
3P,_,, were 5480 A and 5380A respectively. On account of the high internal 
consistency of the calculations, Garstang concluded that the identifications of 
Edlén and Waldmeier were incompatible with these results. 

Now the theory underlying Garstang’s calculations is imperfect in one respect. 
It takes into account only approximately, and in an indirect way, the effects of 
configuration mixing, which are known to be large in the configuration 
(1s?) (2s?) (2p). In the present paper I have redone the calculation for Ca xv, 
using an improved theory that explicitly takes account of the major effects of con- 
figuration mixing. The results tend to support the identifications of Edlén and 
Waldmeier. 

The chief difficulty one ordinarily encounters in estimating the effects of 
configuration mixing stems from the multiplicity of perturbing configurations. 
In Section 2 we shall derive some elementary properties of term spectra in 
isoelectronic sequences. We shall find that in highly ionized atoms the problem 
of configuration mixing simplifies considerably. ‘Thus, in Ca xv the only per- 
turbing configuration that needs to be taken into account is (15) (2p). 


2. Theoretical considerations 


2.1. Term energies and complexes.—The term energies of an N-electron atom 
of nuclear charge Z are defined as the eigenvalues of the Hamiltonian 
N h2 Ze 
H(N,Z)=2 (- — V?- =) 
i=1 2m 


v; 


2 
ee. (2.1) 


ivi” 


We shall find it convenient to take as the units of energy and length 

_ mZ*e ne 

— = mz" (2.2) 
N 2 2 

H(N,Z)= = (-ve- =) 42: gt. (2.3) 


i=1 i iis 


In these units 


Ordinarily one treats Z as a numerical constant. In the present discussion, 
however, we shall assume that all radial wave functions and matrix elements are 
explicit functions of Z. This will enable us to survey certain regularities in the 
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structure of isoelectronic sequences that are not apparent in the spectra of 
individual atoms. 

For a fixed value of N the eigenvalues of H(N, Z) depend on Z in a definite 
way. We shall now determine the nature of this dependence. 

When Z-!=0 the right-hand side of (2.3) reduces to a sum of one-electron 
Hamiltonians. The eigenfunctions are antisymmetrized products of hydrogenic 
wave functions and certain linear combinations of these products. Let the symbol 
(n) denote the set of principal quantum numbers (n,, ,, ..., my), taken in any 
order. Then any linear combination of antisymmetrized products, each of which 
has the same set of principal quantum numbers (), is an eigenfunction of (2.3) 
with Z-1=o belonging to the eigenvalue 

E®=— = n;?. (2.4) 
nj in (n) 

From these well-known facts it follows that the eigenfunctions of H(N, Z) 

and the corresponding eigenvalues can be written as power series in Z~?: 


2 
E= (Sr) {EO + EMZ-14 FOZ-2+ ...}, (2.5) 
Ye=%O+Z1PO+ .... (2.6) 


E® in (2.5) is given by (2.4), and ¥%, in (2.6) is a linear combination of 
antisymmetrized products of hydrogenic wave functions, each product being 
characterized by the same set of principal quantum numbers (7). 

To secure faster convergence, one should in practice replace Z in equation 
(2.5) by Z—Z,, where Z, is the shielding constant given by Slater’s rules (8). 

The expansions (2.5), (2.6) are valid for sufficiently small values of Z-1. We 
shall assume that they hold for all Z-4< N-1. 

We shall call a group of terms belonging to a particular set (m) and a particular 
parity pa complex. According to (2.4) and (2.5) the terms of a given complex are 
“intertwined”’: the interval between terms of the same complex is asymptotically 
proportional to Z, while the interval between terms that belong to complexes 
characterized by different sets (”) is asymptotically proportional to Z?. 

We may conveniently designate a complex by the symbol (172°3"...), 
adding the superscript o if the parity is odd. 

A complex contains one or more entire configurations. For example, the 
complex (1723) contains the configurations (1s”) (2s3p), (1s") (293d), and 
(182) (2p38). 

2.2. Configuration mixing in isoelectronic sequences.—F or the matrix representa- 
tion of H(N, Z) we adopt, as is usually done, an LS-coupling scheme. But in 
addition we shall assume that the radial wave functions contain Z explicitly, and 
are so constructed that the diagonal matrix elements have the proper quadratic 
dependence on Z, as given by (2.4) and (2.5).* 

H(N, Z) is diagonal in L, S, and p. Off-diagonal elements connect terms 
belonging to different configurations, as well as terms of different seniority within 
the same configuration. All the off-diagonal terms are asymptotically proportional 
to Z. Using this fact (which follows from the form of (2.3) and the requirement 
that the diagonal matrix elements have the proper quadratic dependence on Z), 


* A method for constructing such radial wave functions has been developed by the writer and 
applied to the complex (1°) (2%). Only an abstract has been published (10). 
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and the usual formulae of perturbation theory, one can easily demonstrate the 
following theorem. The coefficient E® in the expansion (2.5) of a given term 
energy is completely determined by the submatrix whose elements connect 
terms belonging to the same complex. Moreover, in the expansion 


PY nyraus ow ( - Qn)PBLS; (n't)a’ LS © (n't yx’ LS (2.7) 
nla’ 


of an eigenfunction in terms of the basic functions, the coefficients have the 
following asymptotic dependence on Z: 


ioe: dail 

const. Z . (n’)A(n) (2.8) 

const. if (n’)=(n). 

Thus, the coupling between terms of the same kind that belong tothesame complex 

approaches constancy with increasing Z, while the coupling between the terms 

of the same kind that belong to different complexes tends to zero with increasing Z. 
Since a complex contains more than one configuration in general, the relations 


between term intervals predicted by the standard theory will not ordinarily hold 
true even in the limit Zo. 


Q(n)PBLS; (n'l’)x’ LS { 


2.3. Extrapolation along isoelectronic sequences.—Let Ej,ypy denote the energy 
of the term I’ of the complex (n)”._ We define the difference operator A, in the 


usual way : 
Af(Z)=f(Z + 1) -f(Z). 
From equation (2.5) we obtain the equations 
A(Ewypy — Eq’) = Efior — Efe’ + O(Z-*), (2.9) 
and A2(Eqyr — Eqyp'r’) = 2(Efiyer — Efe’) + O(Z-). (2.10) 
Because of the rapidity with which the differences (2.9) and (2.10) approach con- 
stancy, these rules enable one to extrapolate term spectra with high accuracy 


along isoelectronic sequences. The constant term on the right-hand side of 
(2.10) is given explicitly by (2.4). 


3. Application to Ca xv 


The theoretical energy levels of an atom (N, Z) are the eigenvalues of the 
Hamiltonian H(N, Z)+V(N, Z), where H(N, Z) is given by (2.1) and V(N, Z) 
represents the spin-dependent interactions, which cause the splitting of terms 
into levels. 

In view of the preceding discussion we may reasonably assume that in Ca xv 
the configuration (1s®) (2s*) (2p) is coupled much more strongly with the con- 
figuration (1s*) (2p*), which belongs to the same complex, than with any other 
configuration. On the simplifying assumption of zero coupling between the 
complex (172*) and other complexes, the energy levels of (122*) are the eigenvalues 
of the matrices (3-1 abc). 

The notation is standard* except for the symbol ®, which represents a 
numerical multiple of the Slater integral G,(2s, 2p), the symbol LZ, which repre- 
sents a linear combination of Slater integrals Fy, and the asterisks, which serve 
to distinguish corresponding parameters in the two configurations. 


* The coefficients a and b have not yet been calculated. Their values are unimportant in the 
present calculation, where, for the sake of simplicity, they have been set equal to 2.8 and 1 
respectively. 5 
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2p? 2p* 
“Sp, 1D, “ap, 1D, 

[?s 5f—6-4n 271/(¢-+2-8n) ® ° 
2p? 

| "Dz 2-1/2 +2-8n) 6F, ° 0) (3-1 a) 

p,| ; L—-st*—3-4n* | —2-0C* +09") 
2 a 
° 1D, ° ® —2-1/(f* + an*) L+6F, 

2p? ?P, 2p‘ 3p, 
(3.1 b) 
2p**P, ® L+-5o*+5* 
2p? 2p* 
“ap, 1S, “3p, 1S) 

(*Po —l+8y —2(¢ +7) ® ° 
2p? , 

[3So —21/4(f+m) 15F, ° 20 

(3.1) 

3p, © ° L+{*+2n* 21/(¢* +bn*) 
2 4 
Pp 1S, > 20 21/#(£*-+ bn*) L+15F,* 




















The nine independent intervals between the eigenvalues of (3.1) involve only 
eight undetermined parameters. Unfortunately, the experimental data* for 
the configuration (1s*) (2p*) are so scanty that one cannot determine the values of 
all eight parameters in Ca xv by extrapolation along the isoelectronic sequence. 
However, by a simple theoretical argument we can reduce the number of in- 
dependent parameters from eight to five. 

According to Slater (8), the shielding constant of a 2s electron for a 2p electron 
is the same as that of a 2p electron. A 2p electron in the configuration (1s?) 
(2s?) (2p?) should therefore ‘‘ see’ virtually the same electrostatic field as a 2p 
electron in the configuration (15) (2p*). Hence to a good approximation we may 
a F,*=F,, (*=¢, n*=7. (3.2) 
This approximation should be particularly good in a highly ionized atom like 
Ca xv, since small changes in the shielding will not alter the total electrostatic 
field appreciably. 

We shall use Feenberg’s (9) fifth-order perturbation theory to calculate the 
three smallest eigenvalues of the matrices (3.1). In Ca xv the results are accurate 
to about one part in a thousand. Since the calculation as a whole is accurate 
only to about one part in a hundred, this method of solving the secular equations 
is precise enough for our purpose. 


* All the experimental data used in this paper are taken from Charlotte Moore’s compilation (11). 
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We obtain: 
E(®P2)=0'5¢ — 6-4n 


(¢ + 2-8n)?/2 











Be (+28n)/2 oe? ) 
(°- rap seatreg) — (2st Ltt 
@2 
ef — aan. bt2enyis 
(z-05¢ oR" Teed teas 














({ + 2-8n)?/2 @? ) 
“i (0 ottiaiti 6F,-05¢+6-4n L—l+3y 
20°(¢ + 2-8n)?/2 
* 6F-o 56+ 6-4n)(L—€+3)(L+6F,—0-5£+ 6-47)’ 


©? 
@2 
L+o5¢+5n— (-0'5¢+8— lt) 





(3.3 a) 


(3-3 b) 





E(@P,)= —0'5£ + 8n— 


ECP,)= ~¢+8y a 
2(f+7 
; 40? 2(¢ +n)? 0 
(5%. LiF atm) ¥ (-£+81- 15F y+ 2-87 r¥at=6) 
@2 
2(¢ +n)? 2(¢ +n)? - 
(u+¢+20- L+15F,+0— x) re (-£+8- 15F,+l—8y L+at=) 
2 x 20? x 2(f +n)? ( 
* (5F,+¢—8)(L+al— 6)(L+15F,+f—8n)" ~— 























We shall also need the formulae 
2 





E(2p**P,)=L +0504 54+ @2 ’ (3.4) 

L+o-5¢+5n— (-0'5¢+8)- rt) : 

1 al 

E('S) =15F,— L440’ (3.5 a) 

EP) = — oaek. (3-5 b) 
L+@?/L’ ‘ 


@2 
E(D)=6F,— LyOUL" (3.5 ¢) 


From (3.5) we obtain 
DP=E('D) — EP) =6F,, (3.6) 


@2 4 I 
SD=E('S)— E('D) =9F,— yA (; +4077? Te ni): (3-7) 





By (3.4) and (3.3 a) i 
E(2p**P,) — E(2p?*P,) = 1+. £—4n + ~ oo: (3.8) 
L+l—3n+ 





L+¢—3n 
47 
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Garstang has derived values of SD, PD, ¢ and » for the first ten members of 
the isoelectronic sequence, and, by extrapolation, for Ca xv as well. The extra- 
polation is probably quite accurate, since the approximations underlying the 
calculation are valid for the early members of the sequence. One could probably 
derive improved values of ¢ and » for Si 1x and P x by using the more accurate 
formulae given above; but Garstang seems to have given little weight to these 
ions in extrapolating the values of and. We shall therefore use his extrapolated 
parameters in the present calculation. This step not only saves labour, but also 
puts the essential differences between Garstang’s calculation and the present one 
in the sharpest possible focus. 

Garstang’s values of SD and PD are well represented beyond Ne v by the 


eens SD = 5250 (Z—3.25)—1785 cm", } 
DP = 4660 (Z — 3.25) — 1930 cm”. 
The shielding constant, 3.25, is that given by Slater (8). The smallness of the 
constant terms in these expressions justifies our assumption of zero coupling 
between the complex (12) and other complexes. 
Experimental values for the interval E(2p**P,)— E(2p?*P,) are available for 
O 11, F tv, Ne v and Na vi. Garstang’s calculations give SD, DP, ¢ and » for 
these ions. We can therefore solve equations (3.7) and (3.8) for ® and L. The 
derived values of L are well represented by the formula 
L=63 175 (Z—3°25) —21 380 cm”. (3.10) 
Hence L =1 036 800 cm"! in Ca xv. (3.11 a) 
Garstang’s extrapolated parameters are 
6F,= DP=76 400 cm“, 
SD = 86 300 cm“, 
C=24210 cm}, 


(3-9) 


in Ca xv. (3.11 b) 


n= gicnr}, 
Fro: : (3.11 a, b) and (3.7) we obtain 
@®=101310cm™! in Caxv. (3.11 c) 
Finally, we insert the values (3.11) in the equations (3.3) and obtain 
E(*P,)= — 2639 cm, 
E(®P,) = — 20966 cm“, (3.12) 
EP.) = — 38 667 cm". 
The intervals are 
E(10)=17701 cm™!=5650A (observed: 56944), 
E(21)=18 327 cm™!=5456A (observed: 5445 A). 


\ (3-13) 


In view of the approximations we have made, the agreement is satisfactory; the 
results tend to confirm the identifications of Edlén and Waldmeier. 


The research reported in this paper has been supported by the Geophysics 
Research Directorate, Air Force Cambridge Research Center, under Contract 
AF 19(604)-146 with Harvard University. 

Harvard College Observatory, 


Cambridge, Massachusetts ; 
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ON THE IDENTIFICATION OF Caxv AND Axiv IN THE SOLAR 
CORONA 


Bengt Edlén 
(Communicated by F. J. M. Stratton) 


(Received 1955 February 9) 


Summary 
The interpretation of the coronal wave-lengths A 5694 and A5446 as 
transitions in Ca xv is shown to be compatible with known isoelectronic 
data on level intervals. 'The original suggestion concerning A XIV is revised 
and the identification transferred from A 4359 to the new coronal line A 4412. 





The interpretation (1, 2) of the coronal lines A 5694 and A 5446 as transitions 
within the ground term 2s?2p?3P of Caxv has been subject to much discussion 
with arguments for and against the identifications (e.g. 3, 4). ‘The present 
note is intended to show in a simple and direct way what evidence the known 
data in the isoelectronic sequence actually provide as far as the size of the intervals 
is concerned. 

The data are collected in Table 1, where columns 2, 3 and 4 give all usable 
observed values of the relevant intervals. In some cases the values differ slightly 
from those in the literature (5), which is due to re-examination of existing 
or, in the case of Al vii, utilization of unpublished measurements. ‘The values 
given in brackets for Sirx and Px have been extrapolated in the present 
investigation*, and the values given for Caxv are the wave-numbers of the 
coronal lines. 

TABLE | 


Isoelectronic data for the configuration 2s*2p* 








lon *P.—*P, | *P,-*P> | *D.-*P, X 2-1 X10 Co-1 Ci-0 
cm~! cm! cm~! 

Ca 27°03 16°36 | 10149 | 0°0034 | 0°9986 | 0°5023 27°07 32°57 
Nil 82:2 49° 15184 | 0:0069 | 0-g971 | 0°5046 | - 8274 97°3 
Ol 193°0 113°6 19967 | 0°0122 | 0°9949 | 075081 194'0 223°6 
F Iv 388-2 225'2 24628 | o-0198 | 0°9917 | 0°5131 391°4 438°9 
Ne v 698 414 29182 | 0°0300 | 0°9874 | 0°5198 707 797 
Na VI 1161 696 33620 | 0°0432 | 0°9817 | 0°5284 | 1183 1317 
Mg vII 1812 1127 37982 00600 | 0:9744 | 0°5392 | 1860 2090 
Al vil 2704 1716 ibe 00806 | 0°9654 | 0°5522 | 2801 3107 
Si 1x 3864 (2563) sae 0'1054 | 0°9542 | 0°5677 | 4049 (4515) 
Px (5345)  |(3705) ‘6 0°1350 | 079408 | 075858 |(5681) |(6325) 
Ca xv 18355 17556 _ 0°3673;| 0°82974| 0°71323} 22121 24615 



































* This extrapolation confirms the suggestion by Garstang (6) and Naqvi (7) that the published 
values for P x need revision. 
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The following treatment is, in principle, identical with that used in my 
original paper (1). It is based on the intermediate-coupling formulae which 
express the level intervals as functions of the parameters ¢ and F, (see 8 and 9). 
Introducing the coupling parameter x =(/5F, and the following abbreviations 

a; ~=(?P;—9P,)/f, \ 
B=(*D,—*P2)/5 Fs, 
we obtain from the intermediate-coupling formulae the relations 


ao er ee eee ae 
XMea= +X (2 tps 2x?) ’ 


(1) 


3 3 2 ; 1/2 

yr gan S41 +2 x+x) ; (2) 

36 6 9 i 

=[— — -y+- 2 ‘ 
‘ (; i | 
The parameter x has been defined in the present case by identifying the 
observed ratio R,=(?P,—*P,)/({D,—*P,) with its theoretical expression 
(Ro= x%:_»/B) as a function of x according to formulae (1) and (2).* The 


x-values obtained in this way for the 7 ions from C1 through Mgvii are 
represented within their error limits by the empirical formula 





\/x = 0°04080Z — 0:1022 + 0°048/Z, (3) 
which yields the y-values in col. 5 of Table I and, by way of formulae (2), the 


«’s in cols. 6 and 7. Finally, the ¢-values in the last two columns of the table 
are obtained by dividing the observed intervals by the corresponding «: 


Coy =(*P2—*P))/a24, \ (4) 
Cio = (FP, — 9 Po)/a1-0- 

The difference between (,_, and ¢,, is a measure of the inaccuracy of the 
intermediate-coupling formulae, and is due to their neglecting spin-spin 
interaction and the perturbing influence of other configurations.+ ‘This 
deficiency, however, is of no consequence for the present purpose, as long as 
the parameters obtained are suitable for isoelectronic comparison and extra- 
polation. 

The results are best demonstrated by a graph (Fig. 1) showing the variation 
of ¥/¢ with atomic number. A suitable linear function of Z is subtracted to 
get an adequate scale. As the figure shows, the curve for ¢, _, is definitely pointing 
towards the value which corresponds to 45446 in Caxv, providing a clear 
evidence in favour of the identification. For ¢, 4 the evidence is less definite 
because the points are considerably more affected by observational errors, this 
being due both to larger errors in the intervals themselves and to the smaller 
value of «. The graph shows in any case that an interval in Ca xv corresponding 
to A 5694 is in no conflict with the isoelectronic data as far as they are known. 

The curves in Fig. 1 can be represented by simple interpolation formulae 
which for some purposes may be more useful than a graph. For ¢,_, the values 


* For convenience the «’s, B and y%»-9/B were calculated to 5 decimals for x= 0-00 (0°01) 0-80 
and tabulated. The tables could be duplicated if they should turn out to be of interest to 
others. 

{ I am grateful to Dr Garstang for a discussion on this point, 
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through the whole sequence, including Ca xv, are represented within the limits 
of observational errors by the formula 


WV bg =0°70124Z — 18118 — 0:2889/(Z — 3-48). (5) 


Expressions of the same form, AZ— B—C/(Z—a), were tested in the course of 
a comprehensive investigation (unpublished) and found well adapted to represent 
4/C, F, and 4/ x for all configurations 2s"2p" and 3s"3p". A theoretical 
justification for this form is implicit in a paper by Layzer (10).* 

The conclusion of this discussion is, then, that the identifications of A 5694 
and A 5446 with °P,_, and *P,_, of Caxv are perfectly compatible with observed 
isoelectronic data. ‘Taken together with other evidence (4) this would seem to 
remove further doubts regarding the interpretation of these lines. 
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Kic. 1.—Graph of #“£—0-71(Z—4) referring to the intervals 3P,_, and *P,_4 of the configuration 
2s*2p*. 


It is of interest to mention also the following point which, in a way, provides 
further evidence for the same interpretation and, at the same time, removes the 
last question-mark from the list of coronal identifications. After the recent 
discovery by Lyot and Dollfus (11) of the coronal line A 4412-4 it appears that 
the original, tentative identification of A 4359 with A xIv 2s?2p?P5,.—?Pj). should 
be revised and the assignment transferred to the new line. This fits the iso- 
electronic data definitely better. Moreover, the characteristic appearance of 
A 4412 is similar to that of A 5694 and A 5446, the three lines forming a group 
of their own, precisely as one should expect from the fact that both A xiv and 
Ca xv are exceptional with regard to ionization potentials. 


* 1 am much obliged to Dr Layzer for his kindness in sending me a manuscript copy of his 
paper. It incited me to publish the present note. 


Department of Physics, 
University of Lund, 
Lund (Sweden) : 

1955 February 5. 
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